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ABSTRACT 


The  experiments  described  in  this  report  were  developed  to  test  and 
evaluate  prototype  tools  and  advanced  techniques  to  predict  geologic  con¬ 
ditions  on  the  sea  floor.  Emphasis  has  been  placed  on  determining  the 
mass  physical  properties  of  sea-floor  sed iments,  and  rocks  using  geophysi¬ 
cal  tools  for  the  identification  of  associated  parameters.  A  direct  cur¬ 
rent  resistivity  system,  fashioned  after  the  land-oriented  Wenner  array, 
was  successfully  developed.  It  tvs  been  shown,  through  theoretical  stud¬ 
ies,  that  it  is  feasible  to  use  electromagnetic  sounding  to  measure  conduc¬ 
tivity  and  thereby  deduce  porosity  at  shallow  depths  below  the  sea  floor- 
water  interface.  A  seismic  shear  wave  generator  prototype  tool  was  built 
and  successfully  tested  to  produce  the  transverse  wave  form  in  saturated 
marine  sediments.  The  reflectivity  experiment  to  classify  sea-floor  sedi¬ 
ments,  by  a  ship  underway,  has  verified  previous  field  measurements  by 
Breslau  (1967)  and  others.  A  relatively  high  degree  of  correlation  links 
acousitc  absorption  (bottom  loss)  to  parameters  of  porosity,  density, 
median  grain  size,  fine-grained  sediment  types  (silt  and  clay)  and  sorting 
coefficients.  These  parameters  are  shown  to  have  a  decided  interdepen¬ 
dence,  and  more  importantly,  provide  necessary  information  for  practical 
application  to  engineering  problems  relating  to  underwater  operations 
and/or  excavation  of  the  sea  floor. 
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I.  INTRODUCTION 


This  report  is  a  compilation  of  investigations  performed  by  four 
individuals  working  in  two  separate  fields  of  endeavor.  The  overali 
project  objective  is  to  determine  the  mass  physical  properties  of  sea¬ 
floor  sediments  and  rock  types  to  predict  the  geological  implications 
involved  in  excavation  of  the  sea  floor  and  subsurface  tunneling. 
Studies  are  d  Irected  to  defining  the  geophysical  methods  most  appllcabl 
to  the  problem.  The  two  fields  of  investigation  encompass  application 
of  electrical  and  acoustical  methods,  and  determining  whether  the  tools 
and  techniques  representing  these  methods  have  application  In  a  marine 
environment. 

The  In  situ  measurement  of  sea-floor  properties,  using  electrical 
prospecting  methods,  has  shown  promise  for  determining  varying  conduc¬ 
tivity  within  the  near  surface  geologic  section.  The  work  described  by 
J.  H.  Beyer  and  R.  Corwin  Indicates  that  conductivity,  as  a  function  of 
depth,  may  be  determined  by  direct  current  resistivity  measurements  and 
according  to  Beyer,  It  Is  theoretically  feasible  to  employ  electro¬ 
magnetic  sounding  to  determine  sediment  conductivity.  Corwin  has 
further  shown  that  land-oriented  direct  current  resistivity  tools  may 
be  modified  for  use  In  the  marine  environment. 

Acoustic  experiments  In  this  report  were  designed  and  tested  by 
T.  Hildenbrand  and  B.  Barnes,  HI Idenbrand's  investigations  included 
evaluation  of  mathematical  models  to  ascertain  whether  shear  waves 
could  be  generated  and  theoretically  detected  In  saturated  material. 

He  developed  a  prototype  Instrument  for  generating  shear  waves  to  test 


His  theory.  The  results  from  field  tests  in  marine  environment  confirm 
the  presence  of  the  transverse  wave  form. 

Based  on  previous  investigations,  it  has  been  proven  that  acoustic 
echoes  reflected  from  the  sea  floor  and  substrata  are  Influenced  by  the 
geologic  nature  of  the  media  and,  furthermore,  It  is  reasonable  to 
assume  that  the  respective  attenuation  and  energy  loss  may  be  correlated 
to  parameters  which  depict  the  mass  physical  properties  of  sea-fioor 
sediments.  The  work  by  Barnes  consisted  of  a  field  exercise  to  verify 
this  latter  assumption  and  to  apply  the  results  to  a  practical  engineer¬ 
ing  problem. 

In  each  task  of  the  overall  project  a  systematic  research  approach 
was  applied;  first,  with  an  Inventory  of  previous  work  in  the  area 
followed  by  an  analysis  of  the  theories  applied  to  a  marine  situation 
and  last,  where  applicable,  prototype  development  to  test  the  feasibility 
for  application.  Much  effort  was  made  to  personally  consult  prominent 
authorities  with  years  of  expedience  in  each  of  the  respective  fields. 
These  individuals  are  appropriately  referenced  In  the  text  and  the  list 
of  references. 

The  philosophy  for  first  approaching  the  sea  floor  interface  before 
considering  the  prediction  of  subbottom  properties  has  resulted  from  a 
preliminary  inventory  of  previous  work  accomplished  to  measure  the 
properties  of  the  sea  floor  and  a  knowledge  of  its  practical  application, 
it  was  discovered  very  early  how  little  has  actually  been  accomplished 
toward  development  of  tools  and  techniques  for  making  measurements  of 
sea-fioor  sediments  and  rock.  Several  "firsts"  In  this  regard  have 
resulted  from  the  effort  to  date.  These  are  reported  in  detail  in  the 
text . 
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As  research  and  development  progresses,  the  objectives  wili  be 
extended  to  include  the  lithologic  variations  of  both  consolidated  and 
unconsolidated  materials  at  depth  as  well  as  on  the  sea  floor.  The 
ability  to  predict  the  geologic  situation  in  advance  of  borehlil  drilling 
and/or  excavation  of  the  ocean  bottom  and  beneath  the  sea  floor  is  the  ulti¬ 
mate  objective  of  each  of  the  tasks  described  in  this  report. 


2.  DETERMINATION  OF  SEOIMENT  PROPERTIES  BY  ELECTRICAL  METHODS 

2.1  Relationship  Between  Electrical  Resistivity  and 
Other  Properties  of  Marine  Sediments 

Electrical  conduction  In  saturated  unconsolidated  sediments  occurs 
through  the  Interstltl tal  water,  the  mineral  grains  being  several  orders 
of  magnitude  more  resistive.  The  resistivity  Is  determined,  then,  mainly 
by  the  amount  of  water  present,  Its  salinity,  and  by  the  way  the  water  is 
distributed  through  the  sediment.  The  volume  percentage  of  pore  space, 
or  porosity  <f>  of  the  mater  it i  determines  the  maximum  amount  of  water 
which  can  be  present.  Thus,  In  a  saturated  sediment  the  resistivity  Is 
primarily  a  function  of  the  salinity  of  the  interstitial  water.  In  order 
to  compare  the  resistivities  of  various  samples  It  Is  necessary  to 
normalize  the  values.  This  generally  is  done  by  calculating  the  formation 
factor,  F  (Archie,  1942),  where 

p 

P  =  bulk  resistivity  of  sample _  „  s 

resistivity  of  interstitial  water  Pw 

Kul lenberg  (1952)  showed  that  for  some  deep-sea  sediments  the  difference 
in  salinity  between  the  interstitial  water  in  the  sediments  and  the  bottom 
water  Is,  on  the  whole,  rather  small,  so  that  Pw  Is  easily  obtainable. 

The  formation  factor  is  a  complex  function  of  grain-size  distribution, 
grain  ihape,  packing  arrangement  (see  Section  2.3.3).  and  the  degree  of 
cement  )tion;  or  looked  at  in  another  sense,  formation  factor  is  a 
function  of  porosity,  permeability,  and  capillarity  (Schopper,  1966), 
all  of  which  determine  the  manner  by  which  electrical  current  passes  through 
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the  material.  in  unconsolidated  sediments  the  pore  spaces  are  well 
interconnected  so  that  the  formation  factor  is  primarily  dependent 
upon  porosity. 

An  empirical  relationship  between  formation  factor  and  porosity 
has  been  developed  by  Archie  (1942): 

F  -  #_m 

where  he  believed  m  to  be  a  function  of  cementation.  However,  Erchui 
and  Naccl  (1970  have  demonstrated  that  for  clean  laboratory  samples  m 
ranges  from  i.6  for  sand  to  3.28  for  sodium  montmori i ionite,  in  measur 
Ing  unconsolidated  marine  sediments  Kermabon  et  ai.  (1969)  have  found 
a  more  limited  range; 

m  =  1.8  to  2.0. 

Wyille  and  Gregory  (1953)  demonstrated  that  Archie's  Law  should 
be  modified  to 

_  ,  "  m 

F  -  a<P 

where  a,  an  empirical  constant  In  the  range  0.4  to  i.4,  Is  a  measure  of 
the  effective  path  length  through  the  material,  in  unconsolidated 
marine  sediments  the  constant  a  is  near  unity. 

Thus  It  appears  that  a  knowledge  of  the  formation  factor  of  an 
unconsolidated  marine  sediment  couid  yield  a  reliable  estimate  of 
oorosity.  From  2,500  measurements  made  on  ocean  bottom  cores,  Kermabon 
et  ai.  (1969)  dispiay  data  Indicating  that  the  formation  factor  is 
sufficient  to  predict  porosity  to  +  6  percent.  Porosity,  or  alterna¬ 
tively  void  ratio,  e,  where 


is  quite  useful  for  estimating  engineering  properties  of  sediments.  For 
example,  in  some  experiments  with  clay  samples  Vees  and  Wlnterkorn  (1967) 
demonstrated  an  inverse  straight  line  relationship  between  the  logarithm 
of  shear  strength  and  the  void  ratio. 

On  the  basis  of  empirical  analyses  of  unconsolidated  marine 
sediments  it  Is  apparent  that  porosity  can  be  well  correlated  with 
other  pronertles  of  the  material.  Hamilton  (e.g.,  1970a,  1970b)  has 
derived  regression  equations  which  describe  low  scatter  plots  of 
porosity  versus  bulk  density,  acoustic  impedance,  the  product  of  density 
times  velocity  squared,  Rayleigh  Reflection  Coefficient,  and  bottom  loss. 
Hamilton's  (1970a)  data  also  demonstrate  a  loose  correlation  between 
porosity  and  mean  grain  diameter. 

2.2  The  Effect  of  Sulfate  Ion  Removal  on  the 
Conductivity  of  Seawater 

As  has  been  discussed  above,  the  resistivity  of  unconsolidated 
marine  sediments  is  very  highly  dependent  upon  the  resistivity  of  the 
interstitial  water,  if  resistivity  measurements  are  to  be  made  for  the 
purpose  of  determining  the  porosity  of  the  sediments  It  Is  necessary  to 
know  whether  or  not  the  resistivity  of  the  Interstitial  water  Is  liable 
to  deviate  significantly  from  that  of  the  overlying  bottom  water.  One 
facet  of  this  problem  deals  with  compositional  changes  of  the  Inter¬ 
stitial  water  with  time.  As  a  result  of  interaction  of  the  Interstitial 
water  with  the  mineral  grains  and  plant  and  animal  life  in  the  sediment, 
ions  will  be  selectively  added  to  or  removed  from  the  wat&r.  This  wlil 
alter  the  ionic  concentrations  and  therefore  the  resistivity  of  the 
interstitial  water. 
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An  experiment  has  been  conducted  to  determine  the  effect  upon 
tne  electrical  conductivity  (the  inverse  of  resistivity)  of  seawater  of 
the  removal  of  the  sulfate  ion.  Of  the  Ions  in  standard  seawater, 
sulfate  has  the  fourth  highest  concentration,  following  chlorine, 
sodium,  and  magnesium. 

Artificial  seawater  was  made  according  to  the  recipe  of  Lyman  and 
Fleming  (1940),  and  was  found  to  have  a  salinity  of  34. 87  parts  per 
thousand  using  a  Hytoch  Portable  Laboratory  Sal i nometer  to  measure 
conductivity,  and  the  international  Oceanographic  Tables  ( 1 966)  to 
calculate  salinity  from  conductivity.  An  aqueous  solution  of  barium 

chloride  was  used  to  precipitate  the  sulfate  Ion  from  the  sea  water  as 
barium  sulfate: 

SO4  +  Ba  +2  Ci" — ►BaS0^|  +  2  Ci“ 

The  BaCl2  solution  was  made  to  have  a  Ci"  concentration  the  same  as 
that  of  the  artificial  seawater  so  that  with  the  addition  of  the  BaCi2 
solution,  the  Cl”  concentration  of  the  seawater  remained  unchanged. 
Assuming  that  the  seawater  was  Oittmarian,  the  Cl"  concentration  was 
found  from  Wilde  (1969,  p.  62)  to  be  0.544  formal.  The  Bad2  solution 
was,  then,  0.272  formal  in  Ba++  concentration,  or  0.284  moial,  based 
upon  an  experimentally  determined  density  of  1.045  g/ctr? . 

The  SO4  concentration  In  the  seawater  was  0,0281  formal  (Wilde, 
1969,  p.  63),  which  Is  0.0288  molal  for  a  density  of  1.024  g/cn ?.  Thus, 
to  precipitate  aii  the  SO4  *  from  the  seawater  as  BaS04  It  was  necessary 
to  titrate  the  seawater  with  BaCij  solution  at  a  ratio  of 


Q«0288  moie^ liter  of  seawater 
0.284  mole/i i ter  of  BaCl2  soin 


0.1014  liter  of  BaCi?  spin, 
liter  of  seawater 
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Individual  samples  of  the  seawater  were  titrated  with  different 
quantities  of  the  BaCij  solution  such  that,  respectively,  20,  40,  60, 

80,  and  100  percent  of  the  SO4  =  would  be  precipitated.  The  treated 
seawater  samples  were  then  filtered  and  the  conductivity  was  measured. 
The  data  are  shown  in  Table  2-1.  The  conductivity  ratios  were 
reproducible  to  within  0.3  percent  for  ail  samples. 

The  addition  of  the  BaCl2  solution  to  the  seawater  has  two  distinct 
effects:  it  significantly  reduces  the  SO4  *  concentration  by  precipita¬ 
tion,  and,  to  a  much  iesser  degree,  It  reduces  the  concentrations  of 
all  the  other  ions,  except  Cl”,  by  dilution. 


Table  2- I. — Conductivity  of  seawater  with  SO4  removed. 


Percent  SO4  “ 
removed 

Conductivity 

ratio 

Conductivity 

(mho/cm) 

0 

0.987 

0.0472 

20 

0.966 

0.0462 

40 

0.944 

0.0451 

60 

0.955 

0.0456 

80 

0.926 

0.0443 

100 

0.939 

0.0449 
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The  effects  of  dilution  of  these  ions  can  be  seen  by  considering 
the  equivalent  conductance,  A  ,  where 

A  =  - ^ - - 

equivale.'its/cnv 

and  a  is  the  electrical  conductivity.  As  Is  described  by  Tobias  (1970), 
the  equivalent  conductance  of  an  electrolyte  at  Infinite  dilution, 

Ao  ,  depends  upon  the  contributions  from  each  of  the  ionic  species 
present,  so 

Ao  ■  A  o  +  A  o 

where  Ao  and  Ao  are,  respectively,  the  summed  conductances  of 
the  cations  and  anions  at  Infinite  dilution.  Seawater  Is  a  complex 
concentrated  solution,  but  clearly,  as  the  Ionic  concentrations  are 
reduced  by  dilution  the  conductivity  will  decrease.  This  Is  demon¬ 
strated  by  the  data  in  Table  2-1. 

The  Increase  In  conductivity  which  was  observed  in  proceeding  from 
80  to  100  percent  removal  of  SO^-  could  result  from  the  addition  of 
BaCi 2  solution  beyond  the  point  of  total  S04=  removal.  In  this 
eventuality  Ba++  Ions  ,.ould  remain  In  solution,  with  the  tendency  to 
increase  the  conductivity. 

The  anomalous  Increase  in  conductivity  at  60  percent  removal  of 
S0^  Is  difficult  to  Interpret.  Assuming  that  It  Is  real,  and  not 
the  result  of  an  error  in  th*2  laboratory,  the  explanation  must  be  that 
the  ratios  of  the  complexes  of  SO4  in  seawater  vary  peculiarly  In 
this  concentration. 
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Based  upon  an  extrapolation  of  work  done  by  Longsworth  (1932) 

on  the  cation  transference  number  —  of  KC I ,  Park  (1964)  calculated 

the  partial  equivalent  conductance,  7T  ,  of  major  ions  in  seawater. 

The  results  are  shown  in  Table  2-2.  (Park1 def ined  equivalent  conduc- 

2  2 

tance  in  terms  of  cm  /ohm  rather  than  cm  /ohm-equivalent.  A  thorough 
2/ 

Table  2-2.—  --  Partial  equivalent  ionic  conductance  of  major  ions 
in  seawater  at  23°C,  assuming  the  cation  transference  number  of 
KCI  in  seawater  is  0.49. 


Cation 

A  + 

Anion 

I- 

+ 

Na 

31 

Cl" 

59 

K+ 

57 

II 

cf 

to 

21 

1  u-4+ 

2  Mg 

i3 

hco3‘ 

15 

Z  Ca4-*- 

19 

Br" 

63 

i  Sr++ 

21 

II 

ro 

O 

O 

-8 

discussion  of  the  significance  of  this  detail  is  given  by  Connors  and 
Park,  1967).  Park  explained  the  negative  conductance  value  for  C0^_ 
as  resulting  from  complexing  of  this  ion  with  some  cations  in  seawater, 
and  he  obtained  his  values  for  the  degree  of  complexing  from  Garrels 
and  Thompson  ( 1 962 ) .  Yet  Park  did  not  discuss  the  complexing  of  S0^-, 
also  described  by  Garrels  and  Thompson.  Table  2-3  displays  an  analysis 


J  The  transference  number,  t,  of  an  ion  Is  the  fraction  of  current 

carried  by  that  ion.  _.  Z!  u: 

c  j  1 1  < 

*'  =  J  cj  zj  uj 

where  c  is  the  concentration,  z  is  the  valence,  and  u  is  the  mobility. 


2/  From  Park  (1964,  p.  734), 
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ofS04  similar  to  that  performed  by  Park  for  C03=,  and  demonstrates 
that  as  a  result  of  complexing,  the  effect  of  S04=  on  the  conductance 
of  seawater  should  be  almost  negligible.  However,  this  is  not  in 

accord  with  the  value  Park  computed  for  Table  2-2  and  has  apparently 
been  disregarded. 


Table  2-3.-  Calculation  of  partial  equivalent  ionic  conductance  of 

S04  ion  in  seawater. 


Dissolved 

species 

V 

Compos i t ion^ 
(Percent ) 

Equivalent 
conductance 
(  ^  “1  citt) 

Contribution  to 
overall  conductance 

(  ft  ~  1  CITT ) 

so4= 

54 

21 

+  11. 

CaS04 

3 

19  for  j  Ca 

-  0.5 

MgS04 

21.5 

13  for  £  Mg 

-  3. 

Na2S04 

21 

31  for  }  Na 

-  6.5 

Total  +  1. 

Park  continued  by  calculating  the 

contribution 

of  major  ions  to 

the  overal 1 

conductivity  of  seawater. 

This  information  is  shown  in 

Table  2-4. 

3/  From  Garrels  and  Thompson  (1962). 
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Table  2-4. 

--Percentage  contribution 
Its  overal 1  electrical 

of  major  Ions 
conductance 

In  seawater 
at  35%.and  2] 

Ion 

Equivalent/I  Iter*^ 

Welqht  (%) 

Conduc ti 

Na+ 

0.483 

30.7 

wviiuuu  L O 

29. 

.  ++ 

Mg 

0.109 

3.7 

2.7 

Ca++ 

0.021 

1.2 

0.77 

K+ 

0.010 

1.1 

l.l 

r.  +  + 

br 

0.0001 

0.0 

0.004 

Cl" 

0.558 

55.2 

64. 

sof 

0.057 

7.7 

2.3 

HCO3* 

0.002 

0.1 

0.06 

Br“ 

0.001 

0.2 

0.12 

If,  In  spite  of  the  discrepancy  mentioned  above,  these  values  are 
assumed  to  be  correct,  they  can  be  used  to  calculate  the  decrease  In 

conductivity  of  the  seawater  which  would  be  expected  from  the  addition 
of  the  BaClj  solution. 


57  From  Park  O'M,  p.  735). 
5/  From  Harvey  (i960,  p.  4). 
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For  100  percent  precipitation  of  SOt,’,  the  contribution  of  this  ion 
hecon.es  aero,  the  other  ions,  except  for  untier9„  a  dilution  o,  ,0 

percent.  These  ions,  which  account  for  33.7  percent  of  the  conductivity 
of  seawater,  have  their  combined  effect  reduced  to  30.3  percent.  This 
means  that  for  the  unchanged  value  for  Cl'.'  6b  percent,  the  conductivity 
of  the  seawater  should  be  reduced  to  94.3  percent  of  Its  initial  value 
upon  the  addition  of  the  8aC 1 2  solution.  „„„  all  SO,,’  removed,  the 
treated  seawater  should  have  a  conductivity  of  0.0445  mho/cm. 

Computations  similar  to  this  for  lesser  quantities  of  BaCI2 
solution  added  to  seawater  would  yield  conductivities  along  the  dashed 

in  figure  2-1.  This  analysis  agrees  very  „lth  ,„e  experimental 

data. 

It  must  be  noted,  however,  that  at  least  60  percent  of  the  reduction 
in  conductivity  is  due  to  the  dilution  of  the  seawater  sample,  while  less 
than  40  percent  Is  due  to  ,h.  remova,  of  ,h.  sulfat.  The  remove, 

of  the  sulfate  ion  alone  would  reduce  the  conductivity  from  0.0472  to 
0.0463  mho/cm,  a  change  of  2  percent,  which  Is  In  agreement  w„h  the 
theoretically  obtained  value  given  by  Park  (1964)  in  table  2-4.  A 
Chang.  In  conductivity  of  2  percent  Is  well  „i,h,„  r,„ge  ,rror 
Of  virtually  all  electromagnetic  and  resistivity  systems,  including 
resistivity  probes,  and  Is  therefor.  Insignificant.  Thus.  I,  is 
apparent  from  table  2-4,  that  Ionic  replacement  of  only  chloride  and 
sodium  Is  liable  to  result  in  substantive  changes  of  marine  sediment 
pore  water  resistivity.  This  means  tha,  a  probe  designed  to  selectively 
measure  the  concentrations  of  only  chloride  and  sodium  Ions  would  yield 

sufficiently  accurate  Information  for  the  calculation  of  pore  water 
resistivity. 
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(  mho  /  cm  } 


2.3  Resistivity  Heter  for  Seawater  and 
Unconsolidated  Sediment  Samples 

2.3.1  Description 

A  four  electrode  device  has  been  constructed  to  measure  the 
resistivity  of  seawater  and  unconsolidated  saturated  sediment  samples. 
The  sample  Is  poured  into  a  tube  one  inch  in  diameter  and  about  one 
foot  long  containing  the  four  equlspaced  electrodes.  An  electric 
current  Is  forced  through  the  sampl  between  th»  two  electrodes  near 
the  ends  of  the  tube,  and  a  resultant  potential  difference  is  measured 
between  the  two  intermediate  electrodes.  The  resistivity,  p  ,  of 
the  sample  is 

P  "  K  -j- 

where 

K  -  geometrical  constant  of  the  sample  holder, 

I  =  current  Introduced  through  the  current  electrodes, 

E  *  voltage  drop  measured  across  the  potential  electrodes. 

The  term  K  Is  a  constant  for  the  sample  holder,  and  the  current  lean 
be  kept  constant  regardless  of  sample  resistivity.  Thus,  the  sample 
resistivity  Is  proportional  to  the  voltage  drop  seen  across  the 
potential  electrodes.  The  Instrument  Is  calibrated  (l.e.,  the  term 
K/I  Is  determined)  on  the  basis  of  readings  obtained  using  samples  of 
known  resistivity. 

There  are  several  criteria  which  should  be  met  by  the  current 
generating  and  voltage  detecting  circuits.  The  current  source  should 
meet  the  following  specifications: 
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1.  Square  wave  output. 

2.  Constant  current  regardless  of  sample  resistivity. 

3.  Current  output  in  the  range  of  0.1  -  1.0  miiliampere. 

4.  Output  current  independent  of  supply  voltage  to  system. 

Direct  current  cannot  be  used  because  the  potential  electrodes 

will  eventually  polarize,  making  the  measured  potential  drop  to  that 
existing  between  the  electrodes  and  the  sample.  Any  a.c.  signal  could 
be  employed  but  a  square  wave  has  the  advantages  that  it  is  easy  to 
produce,  and  when  rectified,  yields  a  signal  which  can  easily  be 
measured  at  iow  frequencies  by  a  meter. 

Figure  2-2  shows  the  circuit  for  the  instrument.  The  square  wave 
generator  Is  described  In  the  GE  Transistor  Manual  (1964).  It  is  a 
modification  of  a  multivibrator  circuit  with  the  fiip-fiop  action 
controlled  by  a  unijunction.  The  frequency  of  the  output  can  conven¬ 
iently  be  changed  by  varying  the  resistance  Rj.  With  a  supply  voltage 
of  +  15  volts  the  square  wave  output  oscillates  between  +  10  voits. 

The  current,  if  kept  constant  regardless  of  the  sample  resistivity, 
does  nor  enter  into  the  determination  of  “esistivity  once  the  instrument 
has  beer  calibrated,  and  therefore  does  not  have  to  be  monitored. 
Co.stant  current  is  maintained  by  the  circuit  which  includes  operational 
amplifier  I.  The  current  I  in  the  feedback  loop  of  the  amplifier  is 
determined  by  I  *  ej/Rj,  where  ej  is  the  input  voltage  to  the  ex.nl  if  ier, 
and  Rj  is  noted  in  figure  2-2.  The  amplifier  increases  the  output  volt¬ 
age  e  to  whatever  value  is  necessary  to  maintain  the  current  I*  The 

o 

amplitude  of  this  current,  considering  the  dimensions  of  the  sample 
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holder,  should  be  In  the  range  of  0.1  to  1.0  mill  I ampere;  below  this 
range  the  output  voltage  Is  too  low  for  reliable  measurement  of  a 
seawater  sample  (  <  1.0  millivolt),  and  above  this  average  electrol¬ 

ysis  of  the  water  might  occur,  producing  bubbles  on  the  electrodes. 

The  constant  current  circuit  Is  set  up  to  maintain  a  current  of 
±  0.68  mill  I ampere  through  the  sample. 

To  make  the  current  I  Independent  of  the  supply  voltage,  two 
zener  diodes  are  connected  back-to-back  across  the  output  of  the  square 
wave  generator.  The  amplitude  of  the  output  voltage  of  the  square  wave 
generator  will  vary  with  the  supply  voltage.  For  this  reason,  were  It 
not  for  the  zener  diodes,  the  voltage  e,  seen  by  the  constant  current 
circuit  (and  therefore  the  current  I)  would  be  dependent  upon  the 
supply  voltage.  The  z.ner  diodes  maintain  constant  at  ±  6.8  volts 
for  any  supply  voltage  whose  absolute  value  Is  greeter  then  10  volts. 

Criteria  which  are  important  for  the  detection  circuit  across 
the  potential  electrodes  are:' 

1.  Hjgh  Input  impedance. 

2.  Amplification  of  signal. 

3.  Common  mode  rejection. 

h.  Rectification  of  a.c.  signal  to  d.c.  for  mater  readout. 

It  Is  Important  that  the  Input  Impedance  of  the  potential  measuring 
circuit  be  high  so  that  the  potential  electrodes  will  drew  virtually  no 
current.  This  avoid,  polarisation  of  the  electrode,,  which  c.n  result  In 
potentials  higher  then  the  voltage  drop  .cross  some  , maple,,  n,, 
potential  electrodes  are  connected  to  operational  amplifiers  2  and  3 
(see  figure  2-21  which  ere  Inexpensive  end  heve  Input  Impedance,  of 
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2.5  x  10^  ohms.  At  this  impedance,  frequencies  down  to  10  hz  can  be 
used  without  significant  electrode  polarization  occurring.  Operational 
amplifiers  employing  field  effect  input  transistors  would  allow  lower 
frequencies  to  be  used. 

Operational  amplifier  4  Is  necessary  for  both  amplification  and 
common  mode  rejection;  It  takes  the  difference  of  the  voltages  of  the 
two  potential  electrodes  (after  these  voltages  have  been  amplified  bv 
the  previous  stage).  The  gain,  A,  of  the  potential  measuring  circuit  is 


and  can  be  switched  to  different  values  as  shown  In  the  box  In 
figure  2-2,  depending  upon  sample  resistivity. 

The  amplified  voltage  difference  between  the  potential  electrodes 
Is  the  Input,  to  operational  amplifier  5,  which  contains  In  Its  feedback 
loop,  an  ammoter  In  a  diode  bridge.  The  bridge  circuit  rectifies  the 
square  wave  current  flow  so  that  it  can  be  measured  by  the  meter,  and 
the  reason  for  locating  the  diode  bridge  In  the  feedback  loop  of  an 
operational  amplifier  is  to  force  the  diodes  to  act  Ideally.  That  Is, 
the  operational  amplifier  raises  the  voltage  across  the  brldgr  circuit 
to  a  level  such  that  the  diodes  will  trigger  and  current  will  flow.  As 
a  result,  very  smai  i  volMge  differences  entering  the  common  mode 
rejection  circuit  can  be  detected,  and  furthermore,  the  detection 
circuit  Is  linear.  The  reading  seen  on  the  ammeter  is  proportional 
to  the  difference  in  potential  of  the  potential  electrodes  and  the 
resistivity  of  the  sample. 
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To  calibrate  the  resistivity  meter,  water  samples  with  different 
resistivities  were  prepared  and  measured  on  a  Hytech  Portable  Laboratory 
Sal i nometer.  The  resistivity  meter  proved  to  be  quite  linear,  as 

I 

designed;  all  known  samples  yielded  calibration  curves  within  2  percent 
of  each  other  over  a  range  of  C.2  to  30  ohm-meters.  For  the  sample 
holder  presently  being  used,  ail  scale  readings  must  be  multiplied  by 
0.^95  to  obtain  the  sample  resistivity. 

2.3.2  Use  of  Marine  Sediment  Resistivity  Meter  for  Reslstlvltv- 
Poroslty  Correlations 

The  sample  holder  is  a  one- Inch  diameter  piastlc  tube  containing 
four  equ is paced  electrodes  at  separations  of  two  inches.  A  field  sample 
is  poured  Into  the  vertically  held  tube  for  measurement.  For  laboratory 
experiments  with  regard  to  the  resistivity  of  marine  sediments  as  a 
function  of  porosity,  the  sample  holder  has  been  modified  In  the  follow¬ 
ing  way:  a  bottomless  graduated  cylinder  (also,  one-inch  diameter)  has 
been  glued  to  the  top  as  an  extension.  With  this  addition,  the  porosity 
of  a  sample  In  the  holder  can  be  determined  simultaneously  with  the 
measurement  of  Its  resistivity. 

Saline  water,  which  will  become  the  Interstitial  water  In  the 
sample,  's  poured  Into  the  sample  holder,  and  Its  volume  (that  in  the 
sample  holder  plus  that  extending  Into  the  graduated  cylinder)  Is  noted. 

At  this  point  the  resistivity  of  the  water  can  be  measured.  Then  the 
dry,  clean  (no  salts)  sample  is  slowly  poured  into  the  holder  so  as 
not  to  entrap  air  bubbles.  The  resulting  change  In  water  level  In  the 
graduated  cylinder  is  the  volume  displace  by  the  sample.  The  position 
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m  the  boundarV  between  the  top  of  the  sediment  and  the  clear  water 
above  yields  the  volume  of  saturated  sediment.  The  porosity,^  ,  of  a 
sediment  sample  is 

0  =  volume  of  void  space _ , 

volume  of  saturated  sample 

Thus,  the  porosity  of  the  sample  in  the  holder  is 

<*>  "  1  -  — V°lume  of  water  a"d  sediment-initial  volume  of  water, 

vo 1 ume  of  saturated  sediment  ' 

The  resistivity  can  then  be  measured  without  disturbing  the  sample. 

Tapping  the  sample  holder,  or  shaking  it  with  a  mechanical 
v.brator,  will  tend  to  pack  the  sediment,  thereby  reducing  the  volume 
of  the  saturated  sediment  while  the  total  volume  of  the  sediment  and 
water  column  remains  unchanged.  The  porosity  can  again  be  calculated 
and  a  new  resistivity  measurement  made. 

2.3.3  Experimental  Results:  Analysis  of  Dolomite  and  Silica  Sands 
Using  the  resistivity  instrument  described  in  2.3.1  and  2.3.2, 
resistivity  and  porosity  measurements  were  performed  on  samples  of 
well  sorted  Ottawa (s i I ica)  and  dolomite  sands  of  the  size  distributions 
given  In  table  2-5.  The  data  is  given  in  Appendix  4.1. 


Table  2-5. -Size  distributions  for  sand  samples 


Mesh 

Particle  size 

in  mm 

i  1 1 1ca 
sand 

Max. 

min. 

mean 

35  x  48 

0.420 

0.297 

0.354 

48  x  65 

.297 

.210 

.251 

100  x  150 

.149 

.  105 

.125 
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Table  2-5  (cont'd) 


_  Particle  s ize  In 


Dolomi te 
sand 


4x6 
9  x  10 
14  x  20 
65  x  100 


4.760 


0.210 


0.841 
0. 149 


0.177 


Figures  2-3  and  2-4  show  the  formation  factors  of  the  samples 
Plotted  against  porosity.  As  is  to  be  expected,  as  the  porosity  of  a 
sample  is  reduced,  the  resistivity  increases. 

In  comparing  figures  2-3  and  2-4,  i,  is  apparent  that  the  Ottawa 
sand  data  is  tightly  grouped  for  all  particle  sizes,  while  the  dolomite 
sand,  for  any  given  porosity,  displays  Increasing  resistivity  as  the 
particle  size  Is  reduced.  (This  phenomenon,  might  be  more  easily  viewed 
in  figure  2-6  where  porosity  is  held  constant.)  Figure  2-5  shows 
magnified  photographs  of  the  dolomite  particles.  Though  the  photographic 
reproduction  Is  poor,  i,  is  obvious  that  the  smeiler  particles  are 
considerably  more  angular  than  the  larger  particles,  suggesting  that 
particle  shape  has  an  effect  upon  formation  factor.  This  Is  reasonable 
because  particle  shape  is  a  major  factor  In  defining  the  configuration 
Of  pore  space,  and  this  determines  the  effective  conduction  path  length, 
or  tortuosity  T,  of  the  material.  Tortuosity  Is  related  to  formation 
factor  and  porosity  as  (Wyllie  and  Gregory,  1953) 

T2  =  <#>  x  F 


POROSITY  % 

t*“”  ^  -•  Ur  0^it. 
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The  dolomite  grains  become  Increasingly  angular  as  the  particle 
size  is  reduced,  so  In  packing  the  particles  to  a  particular  porosity, 
the  tortuosity  Increases  with  decreasing  grain  size.  The  data  displayed 
in  figure  2-6  yields  the  values  shown  In  table  2-6  for  tortuosity. 


Table  2-6.  —  Tortuos  1  ty  values 


Mesh 

<t> 

F 

T 

Dolomite  sand 

4x6 

49 

3.01 

12.2 

9  x  10 

49 

3.22 

12.6 

14  x  20 

49 

3.79 

13.3 

SI  I  lea  sand 

35  x  48 

40 

3.56 

12.0 

48  x  65 

40 

3.80 

12.3 

100  x  150 

40 

3.81 

12.4 

The  fact  that  the  silica  sand  yields  tightly  grouped  curves  In 
figures  2-4  and  2-6  and  only  a  slight  variation  in  tortuosity  suggests 
that  ths  different  sizes  of  sand  particles  are  of  nearly  the  same  shape. 

In  Archie's  Law  (F  -  ♦  the  exponent  m  can  be  found  by  plotting 
the  logarithm  of  formation  factor  versus  the  logarithm  of  porosity;  m 
Is  the  slope  of  the  line  defined  by  the  data.  Figures  2-7,  2-8a,  and 
2-8b  show  the  data  points  collected  for  silica  and  dolomite,  respectively. 
The  data  for  the  silica  sand  is  well  grouped  about  a  line  which  Indicates 
a  value  of  m  of  about  1.3,  which  is  the  value  obtained  by  several 
researchers  for  spheres  (e.g.,  Wyllle  and  Gregory,  1953).  This 

suggests  that  the  grains  of  silica  sand  are  well  rounded  for  all  the 
size?  tested. 
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10  2.0  3.0  4  0 

CONDUCTIVITY,  INTERSTITIAL  WATER  (mho/meter) 

Figure  2-6.  Conductivity  of  sediment  va.  oonduotivity  of  interstitial  water. 
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FORMATION  FACTOR  (  P\  /  PSy, ) 


POROSITY  % 
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The .format ion  factor  versus  porosity  data  points  for  dolomite 
are  not,  as  a  whole,  tightly  grouped  except  for  particular  particle 
sizes.  As  is  shown  In  figures  2-8a  and  2-8b,  as  the  particle  size 

I 

decreases  the  values  of  m  are  about  -1.2,  -1.6,  -1.7,  -1.8  for  the 
sizes  tested.  This  Is  consistent  with  the  evidence  shown  above  that 
the  smaller  particles  are  more  angular  and  will  yield  higher  tortuosity 
values. 

These  tests  of  the  resistivity  Instrument  using  sands  were  quite 
encouraging.  The  curves  shown  In  figure  2-6  of  sediment  conductivity 
versus  interstitial  water  conductivity  are  quite  linear,  as  they 
should  be.  The  average  value  of  formation  factor  shown  in  this  figure 
for  40  percent  porosity  Ottawa  sand  is  3-72.  The  vaiue  determined  by 
Erchul  and  Nacci  ( ! 971 )  was  3.83. 

The  curves  of  formation  factor  versus  porosity  are  quite  satisfying 
given  the  additional  knowledge  of  the  shape  of  the  dolomite  grains  as  a 
function  of  size- 

2.4  Initial  Analysis  of  a  Marine  Electromagnetic  Sounding  System 

Probes  are  presently  in  use  (Kermabon,  et  al.,  1969;  Erchui  and 
Nacci,  l 97 i )  which  measure  the  resistivity  of  seafloor  sediments. 
However,  these  probes  disturb  the  sediment  and  sample  only  a  smaii 
volume  of  the  sediment  through  which  they  pass.  One  method  of 
resistivity  surveying  which  would  allow  the  sampling  of  a  large  volume 
of  undisturbed  sediment  wouid  be  electromagnetic  sounding. 


2.^.1  System  Description 

Th«  electromagnetic  exploration  system  envisaged  employs  as  a  source 
a  horizontal  coll  of  wire  (diameter  about  one-half  meter)  to  which  an 
alternating  current  Is  applied.  At  known  horizontal  distances  In 
plane  of  the  source  coll,  pairs  of  small  receiving  colls  are  located 
Which  measure  the  horizontal  and  vertical  components  of  ,h,  magnetic 
field  strength.  It  Is  essential  that  the  proper  orientations  of  the 
receiving  colls  with  respect  to  the  transmitting  coll  be  ma Inta Ined. 

For  this  reason  they  must  be  secured  to  a  rigid  frame,  which  would 
limit  the  maximum  coil  separation  to  about  5  meters;  a  system  even 
this  large  would  be  rather  unwieldy  to  handle  on  shipboard  and  lower 
to  the  seafloor.  The  coll  separation  must  be  made  as  large  as  possible 
to  maintain  a  useful  depth  of  exploration,  In  this  case  the  upper  I 
to  2  meters  of  bottom  material. 

2.4.2  Analysis 

In  a  homogeneous  medium  there  would  ba  no  horizontal  (radial) 
component  of  the  magnetic  field  In  the  plane  of  the  transmitting  coll. 
However,  If  the  coll  Is  at  the  boundary  of  a  more  resistive  half  space, 
such  as  the  seafloor,  the  field  will  be  distorted,  and  ,h,  ratio  of  ,h. 
horizontal  to  vertical  components  of  the  field  will  Increase  from  zero 
with  the  resistivity  ratio,  or  formation  factor.  Thus,  the  horizontal 
component  of  the  field  Is  sensitive  to  bottom  resistivity. 
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The  amplitude  of  the  vertical  component  is,  on  the  other  hand, 
reiatively  insensitive  because  deviations  due  to  the  bottom  material 
are  oniy  a  smali  fraction  of  the  observed  vertical  field.  However,  the 
out-of-phase  component  of  the  vertical  field  is,  In  the  vicinity  of  the 
source,  due  to  the  damping  effect  of  the  conducting  „.adia,  and  Is  there¬ 
fore  reiatively  sensitive  to  bottom  resistivity. 

The  signal  seen  at  the  receiving  coiis  Is  highly  dependent  upon  the 
frequency  of  the  transmitted  signal.  At  low  frequencies  there  Is  very 
little  Induction  (dB/dt)  at  the  detectors,  and  at  very  high  frequencies 
tfe  field  Is  severely  attenuated  as  a  result  of  the  high  conductivity 
of  the  sea  water  and  saturated  sediment.  The  amplitude  of  the  fieid 
varies  as 

-»/T  ) 

where 

R  =  distance  from  the  source 
a)  =  angular  frequency 
a  =  conductivity  =■  (resistivity)  _i 
y  =  permeabii ity 

Coggon  and  Morrison  ( i 970)  have  calculated  the  fields  due  to  a 
vertical  dipole  (horizontal  ioop)  on  a  homogeneous  sea  floor;  for  a 
coll  separation  of  5  meters  the  frequency  range  which  yields  the 
strongest  response  is  from  300  hz  to  300  k  hz.  At  the  upper  end  of 
this  frequency  range,  however,  It  appears  that  it  would  be  quite 
difficult  to  build  a  transmitting  coli  of  suitable  physical  size. 
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Wfth  a  receiver  at  a  fixed  distance  from  the  transmitter  parametric 
sounding  can  be  accomplished  to  determine  the  response  of  the  sediment 
as  a  function  of  frequency.  A  computer  program  written  by  J.  H.  Coggon  6/ 
has  been  used  to  calculate  the  horizontal  magnetic  fieid  intensity  of 
the  sea  floor;  the  curves  and  sea-fioor  models  are  shown  in  figure  2-9. 

It  is  apparent  that  one-iayer  and  two- layer  sea  bottom  models  yield 
overlapping  families  of  curves.  For  this  reason,  using  only  one 
detector  it  would  be  impossible  to  distinguish  between  a  one-layer 
sea  fioor  with,  say,  formation  factor  F  =  TO,  and  a  two-iayer  case  with 
F  =  5  for  the  top  i  meter  and  F  =  100  beneath  that.  Clearly,  a 
particular  curve  represents  no  unique  case. 

To  overcome  this  difficulty,  it  is  necessary  to  conduct  geometric 
sounding:  changing  the  transmitter-receiver  separation.  This  can  be 
easily  achieved  in  a  marine  electromagnetic  system  by  locating  several 
detectors  at  intervals  along  a  frame  extending  from  the  transmitting 
co i  The  de^tor  nearest  to  the  transmitter  would  have  the  shallowest 
depth  of  penetration  and  would  yield  a  resistivity  estimate  of  the  upper¬ 
most  portion  of  sediment.  This  information  couid  be  used  to  aid  in 
interpreting  the  data  from  successively  more  distant  detectors.  Families 
of  curves  such  as  shown  in  figure  2-9  would  have  to  be  computed  for 
interpreting  the  response.  Alternatively,  interpretation  could  be 
based  upon  the  gradient  of  the  fieid  (actually,  the  difference  of  the 
field  strengths  seen  at  two  detectors). 


It  can  be  noted  from  figure  2-9  that  as  the  resistivity  of  the 

bottom  increases,  the  horizontai  magnetic  field  becomes  less  sensitive 
~  Dept‘  J  als  lienee  and  tnglneer ing,  Unlv.  of  Calif.,  Berkeley 
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to  bottom  resistivity  charges.  For  formation  factors  greeter  than  10 
the  sensitivity  becomes  very  poor.  A  similar  effect  is  noticed  for  the 
imaginary  part  of  the  vertical  component  of  the  field. 

2.4,3  Conclusions 

Further  analysis  of  the  marine  electromagnetic  system  needs  to 
be  done  to  detenr!ne  the  utility  of  measuring  gradients  of  the  horizontal 
field,  and  to  evaluate  more  fully  If  there  is  any  advantage  to  be  gained 
by  measuring  the  Imaginary  part  of  the  vertical  component  of  the  field. 

Presently  there  appear  to  be  two  severe  limitations  on  a  marine 
electromagnetic  sounding  system:  (1)  depth  of  exploration  greater  than 
1  to  2  meters  would  necessitate  a  frame  whose  size  would  be  unmanageable, 
yet  the  importance  of  knowing  coll  separation  and  relative  orientation 
would  seem  to  require  rigid  connection;  (2)  the  system  would  be  rather 
insensitive  to  formation  factors  greater  than  19. 

An  electromagnetic  sounding  system  could  not  distinguish  details 
of  layering  to  the  extent  possible  with  a  resistivity  probe,  but  if 
bulk  properties  of  the  upper  portion  of  sediment  were  desired,  an 
electromagnetic  system  could  prove  quite  useful.  It  would  yield  an 
average  resistivity  over  a  large  volume  of  undisturbed  sediment,  and 
would  reveal  gross  layering. 

2.5  Resistivity  Measurement  by  Direct  Current  Methods 
The  direct  current"  (D.C<)  resistivity  methods  have  been  used 
for  many  years  to  obtain  Information  about  the  structure  and  porosity 
of  subsurface  formations  on  land.  Adapted  to  marine  use,  they  could 
provide  a  valuable  geophysical  tool  for  sea-floor  Investigations. 


36 


Th.re  are  three  principal  methods  of  measuring  subsurface 
resistivity: 

'*  lnductJve  electromagnetic  techniques,  in  which  high  frequency 

alternating  electric  fields  mr*  *.«  •  j 

neias  are  used  to  mduce  currents  and  magnetic 

fields  in  the  ground  (an  application  of  this  is  discussed  in  the 
previous  section). 

2.  Well  logging  techniques,  in  which  a  small  probe  is  lowered 
•nto  a  borehole  or  pushed  Into  soft  sediments. 

3.  "Direct  current"  methods,  In  which  direct  or  low  frequency 

alternating  currents  are  Introduced  directly  Into  the  ground,  and 

potentials  are  measured  by  electrodes  Inserted  into  the  surface  of  the 
ground. 


2.5. 1  Wei  1  Logging  Method 

The  well  logging  method,  as  described  by  Dakhnov  (1959)  and  other 
authors, and  sho»,  schematically  In  f„ur.  2-10,  has  been  used  for  sea¬ 
floor  resistivity  measurements  by  Erchul  and  Nacci  (1971),  Kermabon, 
et  al  (1969),  and  Houma,  et  al  (1971).  Their  results  indicate  that 
the  porosity  and  density  „f  marine  sediments  may  be  obtained  from  the 
measured  electrical  resistivity.  (The  relation  between  electrical  and 

other  physical  properties  of  sediments  Is  discussed  more  fully  l„ 
section  2.1). 

The  use  of  a  logging  probe  gives  a  relatively  precise  determination 
Of  resistivity  as  a  function  of  depth,  but  does  offer  some  drawbacks. 
Measurements  cannot  be  made  while  the  ship  Is  underway,  and  no  measure¬ 
ment  can  be  obtained  If  the  probe  is  unabie  to  penetrate  the  bottom. 

6'so,  the  volume  of  bottom  materia,  sampled  by  the  probe  Is  small,  which 
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,ElectroniC' 


Lower  Sediment  Layer  pz 

A ;  Current  Electrodes 
B  ■  Potential  Electrodes 


Current  I  is  fed  into  the  sediments  through  current  electrodes  A,  and  causes 
a  voltage  AV  ta  appear  between  potential  electrodes  B  Sediment  resistivity 
P  is  then  determined  from  the  formula 


I 

where  K  is  o  factor  depending  on  the  geometry  of  the  electrode  arrangement 


Figure  2-10. 


Well  logging  probe  used  to  measure  resistivity 
of  marine  sediments.' 
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may  lead  to  interpretation  errors  if  the  sampled  point  is  not  representa¬ 
tive  of  the  surrounding  area. 

2.5.2  Horizontal  Array  Method 

Description. 

The  use  of  a  horizontal  array  of  electrodes  (figure  2-11)  is  another 
standard  method  of  resistivity  surveying  on  land.  Mooney  and  Wetzel 
(1956),  Van  Nostrand  and  Cook  (1966),  Keller  and  Frischknecht  (1966), 
Grant  and  West  (1965),  and  Parasnls  (1962),  along  with  many  other 
authors,  cover  this  topic  in  detail. 

Briefly,  a  direct  or  low  frequency  alternating  current  I  is 
introduced  into  the  ground  through  two  (or  more)  electrodes,  usually 
in  the  form  of  copper  or  iron  stakes.  The  potential  difference  on  the 
surface  of  the  ground  is  then  measured  between  two  (or  more)  electrodes 
inserted  into  the  ground.  These  may  be  of  the  same  type  as  the  current 
electrodes  or  they  may  be  of  the  "non-polarizing"  variety.  The  voltage 
AV  measured  between  the  potential  electrodes  is  a  function  of  the 
current  I.  and  of  the  resistivity  and  structure  of  the  ground  below. 

By  expanding  the  array,  more  current  Is  forced  to  flow  at  greater 
depths,  and  Information  may  be  obtained  on  the  change  of  resistivity 
with  depth. 

Many  different  array  configurations  are  used  for  land  surveys. 

The  simplest  Is  the  Wenner  array  (figure  2-1 1),  in  which  ail  the 
spacings  a  are  equal.  This  the  only  configuration  used  In  this  study 
to  date  although  others  must  be  examined  before  the  optimum  array  for 
marine  use  may  be  chosen. 
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Figure  2-11. 
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",t*rPr*tati°"  °f  field  '■  done  by  use  of  .  set  of  curves 

relating  the  ratio  4V/1  to  the  spacing  a  (figure  Us„ally, 

the  ratio  is  converted  to  an  "apparent  resistivity",  ,  by  use  of  the 
formula 

pa  ■  2ira  nr 

The  apparent  resistivity  would  be  egual  to  the  actuai  resistivity 
if  readings  were  made  over  an  Infinite,  homogeneous  half-space  of 
resistr-pty  pa.  As  may  be  seen  from  figure  2-12,  pa  changes  with 
depth  in  a  different  manner  for  different  layered  situations.  Use  of 
the  apparent  resistivity  concept  gives  a  better  guaiitative  fee,  for 
the  resistivity-depth  distribution,  although  theAV/I.  ratio  would  serve 
equally  well  for  quantitative  Interpretation.  this  study,  the 

4V/I-  rat'°  'S  USad  aS  tha  “Tdlnate  for  the  interpretation  curves, 
since  the  magnitude  of  the  differences  this  ratio  wl  „  determine 
whether  cases  of  geological  interest  may  be  distinguished. 

Marine  Appi ication. 

The  use  of  an  expanding  horizontal  array  for  sea-fioor  measurements 
(figure  2-13)  has  a  number  of  advantages  compared  with  the  use  of  a 
probe.  The  "depth  of  penetration"  is  limited  only  by  the  size  of  the 
array  and  the  amount  of  current  available,  so  information  may  be 
Obtained  at  considerable  depths  below  the  seafloor.  Also,  measurements 
may  be  mad.  with  the  ship  underway,  at  a  considerable  savings  in  time 
and  expense.  Even  an  array  on  the  bottom  may,  In  many  cases,  be  dragged 
continuously  across  the  sea  floor  as  done  by  HarKe  (,965).  On  the  other 
hand,  a  probe  gives  a  more  precise  determination  of  resistivity  with 
depth,  and  Interpretation  Is  simpler  because  there  is  no  need  to  refer 
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Ftgure  2-13.  Surface  and  sea-floor  Wenner  arrays 
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to  a  set  of  curves  to  obtain  the  res i st i v i ty--depth  relationship  from 
the  observed  data.  Also,  the  curves  in  some  cases  contain  ambiguities 
which  make  interpretation  more  difficult. 

The  probe  and  horizontal  array  methods,  then,  are  complimentary. 
Since  considerable  work  has  been  done  on  the  probe  method  by  the  authors 
cited  above,  it  was  decided  for  this  project  to  investigate  the 
horizontal  array  techniques  for  their  applicability  to  sea-floor 
resistivity  measurements. 

Some  previous  work  has  been  done  in  this  area.  Marke  (1965)  has 
studied  the  single  layer  case  in  detail,  and  has  dragged  a  fixed  array 
over  the  bottom  off  the  coast  of  Cornwall,  England.  Volker  and  Dijkstra 
(1955)  have  measured  the  chloride  content  of  deep  ground  waters  under 
submarine  sediments  in  the  Netherlands  using  a  marine  Wenner  array. 
Terekhin  (1962)  and  Van'yan  (1956)  have  examined  the  theory  of  the 
multilayered  submarine  case  with  current  and  measuring  electrodes  placed 
on  the  sea  floor.  Other  Russian  workers  have  apparently  published  curves 
for  sea-floor  resistivity  arrays,  and  have  carried  out  resistivity 
measurements  at  sea,  but  we  have  not  yet  been  able  to  obtain  copies  of 
their  work. 

In  order  to  determine  the  utility  of  the  horizontal  array  for  the 
measurement  of  sea-floor  resistivity,  sets  of  curves  similar  to  those 
shown  in  figure  2-12  must  be  developed  for  arrays  deployed  over  a 
layered  bottom,  both  at  the  sea  surface  and  on  the  sea  floor.  A 
comparison  of  these  curves  would  show  whether  the  surface  array,  which 
must  "see"  through  a  highly  conductive  layer  of  sea  water  before  it  can 
sense  the  bottom  conductivity,  would  be  useable  for  practical  surveying 
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(it  is  obviously  much  easier  to  tow  an  array  on  the  surface  than  to  drag 
it  over  the  bottom).  The  work  of  Marke  (1965)  indicates  that  an  array 
on  the  surface  is  much  less  efficient  than  one  on  the  bottom,  but  a 

i 

more  detailed  study,  especially  for  layered  bottoms,  must  be  made. 

Next,  It  must  be  determined  whether  it  will  be  possible  to  select 
or  design  instrumentation  capable  of  measuring  the  desired  quantities 
with  the  accuracy  needed  for  practical  application.  finally,  the 
instrumentation  must  be  selected  (or  designed  and  built)  and  tested  in 
the  marine  environment,  and  the  test  results  compared  with  the  parameters 
measured  on  the  samples  taken  from  the  MMTC  test  range  (section  3.2.8). 
Development  and  Interpretation  of  Curves. 

For  an  array  deployed  on  the  sea  surface,  published  sets  of  curves 
such  as  those  by  Mooney  and  Wetzel  (1956)  for  the  Wenner  array,  or 
Rl jkswaterstaat  (1969)  for  the  Schlumberger  arrangement  are  available. 

For  an  array  on  the  sea  floor,  however,  no  general  sets  of  curves  are 
available.  Terekhln  (1962)  has  developed  the  mathematics  needed  to 
solve  the  most  general  case,  that  of  a  point  source  of  current  on  the 
sea  floor  over  an  N  -  layered  bottom,  but  the  form  of  his  solution  is 
not  ideal  for  computer  calculation  of  actual  values  over  a  bottom  with 
one  or  two  layers. 

Since  three  layers  are  usually  the  most  that  can  be  resolved  in 
practice  on  land,  two  layers  will  probably  prove  to  be  the  upper  limit 
of  resolution  for  marine  studies.  Therefore,  the  mathematics  and 
computer  programs  to  calculate  the  potential  distribution  over  a  two¬ 
layered  bottom,  due  to  a  point  source  of  current  on  the  sea  floor, 
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were  developed.  The  mathematical  analysis  and  the  computer  program  are 
presented  In  section  4.2. 


The  program  calculates  V^/I,  the  voltage  on  the  sea  floor  at 
a  distance  r  from  the  source,  divided  by  the  input  current.  The 
value  of  AV/I  for  the  Wenner  array  may  then  be  easily  found  by  use 
of  the  formula 


where  V2r  is  the  voltage  at  a  distance  r  from  the  source.  Finally, 
the  quantity  AV/I  Is  plotted  against  a  ,  giving  curves  such  as 
those  shown  in  figures  2-14  and  2-15. 


In  order  to  be  of  use  In  the  interpretation  of  sediment  layering, 
the  curves  for  various  cases  of  geologic  Interest  as  obtained  from 
field  data  must  be  distinguishable.  As  discussed  below.  It  appears 
at  differences  on  the  order  of  1  millivolt  may  be  resolved  by  the 
field  equipment  which  has  been  constructed.  It  is  apparent  that  resolution 
is  increased  by  increasing  the  Input  current  I ,  as  this  increases  the 
signal  to  noise  ratio.  As  the  current  Is  increased,  more  subtle  geologic 
cases  may  be  distinguished  from  each  other.  Since  the  system  discussed 
below  will  provide  10  amps  of  current,  differences  In  the  value  of 
AV/I  of  about  0.1  mi  1 1 Ivol ts/amp  should  be  distinguishable  from  the 
field  data.  This  value  will  be  used  as  a  preliminary  criterion  as  to 
whether  two  readings  obtained  in  the  field  may  be  resolved:  if  the 
difference  in  the  value  of  AV/I  for  two  geologic  cases,  at  the  same 
separation  a  ,  is  greater  than  0.1  mill ivol ts/amp,  the  curves  of 
AV/I  vs  a  for  the  two  cases  may  be  distinguished  at  that  point. 
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h  =  IQ  rtt 


Figure  2-14.  AV/I  ve.  a  for  a  surface  and  a  bottom  Wenner  array 
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Figures  2-14  and  2-15  are  examples  of  the  types  of  curves  to  be 
evaluated.  Figure  2-14  shows  the  effect  of  varying  the  basement  layer 
resistivity  while  holding  the  geometry  constant  for  a  hypothetical 
submarine  situation.  It  may  be  seen  that  tlhe  difference  In  the 
resistivity  of  the  bottom  layer  may  be  distinguished  for  values  of 
a  greater  than  I  meter  for  en  array  on  the  bottom  and  for  values  of  a 
gi  eater  than  about  6  meters  for  an  array  towed  on  the  surface.  Figure 
2-15  shows  the  effect,  for  a  different  layered  case,  of  varying  the 
depth  to  the  lowest  layer  and  the  basement  resistivity.  The  higher 
basement  resistivities,  25  and  250  ohm-meters  (ft  -  m)  are  difficult 
or  Impossible  to  distinguish  at  any  but  very  large  values  of  a  .while 
the  lower  resistivity  cases  are  distinguishable  for  values  of  a  ranging 
from  2  to  4  meters. 

Equipment. 

Because  commercially  available  resistivity  units  are  designed  to  be 
used  In  materials  of  high  resistivity,  they  operate  at  high  voltages  and 
low  currents.  A  unit  to  be  used  In  the  low  resistivity  marine  environ¬ 
ment,  however,  must  operate  at  high  currents  In  order  to  give  useful 
values  of  the  AV/ I  ratio,  while  the  low  contact  resistances  Involved 
reduce  the  need  for  high  voltage.  Accordingly,  a  unit  which  uses  an 
automotive  or  similar  storage  battery  as  a  current  source,  and  provides 
square  wave  output  £t  currents  up  to  10  amps  and  at  frequencies  ranging 
from  0.2  to  20  cps.  was  designed  and  built  at  MMTC.  The  wiring  diagram 
for  this  unit  Is  shown  In  figure  2-16.  A  high  Impedance  rectifying 
voltmeter  for  reading  AV  was  also  designed  and 
(figure  2-17). 
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Field  Tests. 


Three  preliminary  field  tests  have  been  run.  For  the  first  test  a 
small  array  with  a  fixed  spacing  a  of  20  cm.  was  lowered  to  the 
bottom  off  the  wharf  at  MMTC  (figures  2-18  and  2-19).  For  such  a 
small  electrode  spacing,  the  bottom  may  be  considered  to  consist  of 
only  a  single  layer  of  resistivity  ps  and  the  water  layer  may  be 
considered  to  be  Infinitely  deep.  In  this  case, 


where 


w 


2f. 


T 


ps  Is  the  resistivity  of  the  (single)  sediment  layer. 
Pw  Is  the  resistivity  of  the  water. 


A  detailed  derivation  of  the  equation  Is  9lv,n  In  Appendix  (Sec.  4.2.1,) 
This  formula  may  then  be  used  to  obtain  sets  of  curves  such  as  that 
shown  In  figure  2-20,  for  a  -  20  cm.  In  this  case,  the  measured 
value  of  AV/I  was  140  ml  I  Hvol  ts/amp,  which  gives  a  value  of 
pS/pW  of  1  *3  for  the  measured  value  of  Pw  -  0.308  8  -  m.  jhe 
value  of  AV/I  appeared  to  be  Independent  of  the  frequency  or  current 
used.  The  chart  recorder  was  used  to  monitor  the  form  of  the  potential 
signal.  It  was  noted  that  the  readings  were  subject  to  D.C.  drift, 
although  the  peak-to-peak  values  of  the  signal  were  reasonably  stable, 
rhe  drift  was  probably  due  to  corrosion  of  the  copper  potential  electodes, 
Indicating  that  nonpolarizing  potential  electrodes  are  a  better  choice 
in  the  corrosive  marine  environment. 
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Figure  2-19.  Teat  aite  -  reaiativity  experiment. 
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( m  mho/ cm ) 


For  the  second  test,  an  attempt  was  made  to  use  an  expanding  Wenner 
array  off  the  same  dock,  deployed  first  on  the  surface  and  then  on  the 
bottom  for  each  spacing  a.  Nonpolarizing  silver-silver  chloride 
potential  electrodes  as  described  by  Von  Ani  (l'962)  were  used,  and 
drift  was  considerably  reduced.  Unfortunately, two  factors  caused  the 
readings  obtained  to  be  extremely  erratic.  First,  the  large  pilings 
and  considerable  amount  of  metal  attached  to  the  pilings  undoubtedly 
distorted  the  potential  fields  severely.  Second,  a  very  strong  tidal 
current  made  It  impossible  to  position  the  electrodes  accurately.  For 
these  reasons,  no  quantitative  Interpretation  was  made  of  the  test 
results. 

The  performance  of  a  towed  array  was  observed  for  the  third  field 
test.  The  array  was  In  the  form  of  a  very  rough  Wenner  configuration 
with  a  ~25  feet  ^figure  2-21).  The  current  electrodes  were  copper 
rods  about  1  Inch  In  diameter  by  about  6  inches  In  length,  while  the 
potential  electrodes  were  of  an  improved  silver-silver  chloride 
nonpolarizing  type  (figure  2-22).  No  attempt  was  made  to  obtain 
quantitative  readings  of  bottom  resistivity  since  the  configuration 
was  not  known  with  sufficient  accuracy. 

An  important  result,  however,  was  apparent  from  the  test:  the 
reading  of  AV  between  the  potential  electrodes,  with  no  power  applied 
to  the  current  electrodes,  showed  negligible  drift  and  noise  levels, 
well  under  1  millivolt  (figure  2-23).  When  power  was  applied  to  the 
current  electrodes,  a  typical  signal  looked  like  figure  2-24.  The 
noise  seen  on  the  signal  was  due  to  problems  in  the  current  source, 
which  are  being  corrected,  and  random  motion  of  the  electrodes,  which 
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Figure  2-21.  Towed  array  field  test. 


*rur  la  seal# 


A  -  Coble  Insulotion 

B  -  Plexiqlos  Housing 

C  -  Conductor 

D  -  Epoxy  Resin  Potting 

E  -  Silver- Silver  Chloride  Element 

F  -  Filling  Solution  -  Soturoted  AgCI  and  KC I 

G  -  Porous  Ceramic  Rod 

H  -  KCI  Crystals 

I  -  Filling  Screw 


Figure  2-22.  Non-polarizing  silver -silver  chloride  electrode 


Figure  2-23.  AV  between  potential  eleotrodee ,  1*0. 


Figure  2-24.  AV  between  potential  electrodes,  l-5.J»amps 
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cotild  be  stopped  by  taping  the  four  lines  together.  Once  these 

problems  are  corrected,  the  signal  should  be  easily  readable  to  less  than 
I  millivolt. 

2.5.3  Conclusions  and  Recommendations 

It  appears  that  the  direct  current  resistivity  method  is 
applicable  to  the  study  of  ..-floor  properties.  Preliminary  field  test 
results  Indicate  that  differences  In  the  value  of  4V  0f  1  millivolt 

may  be  r'3 * S°'V*d  and  Nation  of  curves  such  as  figures  2- lb  and  2-15 
Shows  that  with  currents  up  |0  amps,  differences  of  a  millivolt  or 
more  appear  between  the  curves  iV/I  Vs.  a  for  a,  leas,  of 

the  geologic  cases  studied. 

Recommended  future  work  Includes: 

1.  A  study  of  the  different  array  configurations  to  determine 
which  is  the  most  useful  In  the  marine  environment. 

2.  A  thorough  theoretical  comparison  between  surface-  and  bottom- 
towed  arrays. 

3.  Compilation  of  a  large  number  of  Interpretation  curves  for 
cases  of  geologic  Interest. 

“I.  Fabrication  of  a  multiple  electrode  expanding  array  which  may 

be  streamed  at  the  surface,  towed  over  the  bottom,  or  deployed  stationary 
on  the  bottom. 

5.  Field  tests  over  the  MMTC  test  ranges  In  San  Francisco  Bay 

(section  3. 2. 8.1)  and  Monterey  Bay.  and  comparison  of  results  with  those 

Obtained  by  coring,  drilling,  and  the  other  geophysical  methods  described 
In  this  report. 
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3.  DETERMINATION  OF  SEDIMENT  PROPERTIES  BY 
ACOUSTIC  METHODS 

3.1  Refraction  Methods 

3.1.1  Relationship  Between  Seismic  Parameters  and  Other  Soil  Properties 
Measured  values  of  compress ional -wave  velocity,  shear-wave 
velocity,  attenuation  coefficients  of  these  two  waves,  and  bulk  density 
can  be  used  to  determine  several  properties  of  sediments.  Rigidity, 
bulk  modulus,  Lame's  constant,  Poisson's  ratio  and  Young's  modulus  can 
be  calculated  directly  while  porosity,  mean  grain  size,  and  shear 
strength  possibly  can  be  found  through  empirical  relationships  with  the 
measured  physical  properties.  Some  previous  work  along  these  lines  has 
been  conducted  by  Hamilton  et  ai  (1970  ),  Hamilton  (1970b,  1971),  Horn 
et  al  (1968),  Schreiber  (I968),  Shunway  (I960),  and  Sutton  et  al  (1957). 

3.1.2  Experimental  Methods 

The  determination  of  the  shear-wave  velocity  In  the  upper  layers 
of  marine  sediments  is  difficult  because  shear  waves  will  not  propagate 
in  water  and  consequently  are  transmitted  poorly  in  a  saturated,  porous 
medium.  On  the  other  hand,  the  compress lonal  wave  is  transmitted  readily 
‘n  water  and  in  saturated  sediments. 

To  date,  no  in  situ  seismic  method  has  been  found  that  will 
directly  give  the  shear-wave  velocity  in  saturated  marine  sediments.  In 
fact,  in  earlier  studies  by  Wood  (i94i) ,  it  was  assumed  there  was  no 
transmission  of  a  shear  wave  in  marine  sediments.  It  was  later  found 
by  Nafe  anu  Drake  (1963)  and  Hamilton,  et  ai  (1971)  that  almost  aii 
marine  sediments  have  measurable  rigidities  and  transmit  shear  waves. 
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V i scoe 1 astometers  have  been  demonstrated  to  be  a  means  of  obtaining 
the  shear-wave  velocity.  The  method  involves  measuring  the  real  and 
imaginary  parts  of  the  rigidity  from  torsional  waves  propagating  along 
the  length  of  a  round  rod.  However,  the  method  is  very  dependent  on 
the  temperature  properties  of  the  rod,  and  furthermore,  an  in  situ 
v i scoe las tometer  has  not  yet  been  developed. 

Another  method  used  to  determine  the  shear-wave  velocity  is  the 
use  of  Stoneley  waves  (Davies,  1965.  Hamilton,  1970,  et  al).  Stonnley  waves 
are  generalized  Rayleigh  waves  that  propagate  along  the  interface 
between  two  media  such  as  a  water-sediment  Interface.  After  obtaining 
the  Stoneley-wave  velocity  from  a  seismogram,  the  shear-wave  velocity 
is  found  from  equations  relating  the  shear-wave  velocity  to  the  Stoneley- 
wave  velocity  (Strlck  and  Glnzberg,  1956).  The  disadvantage  of  this 
method  is  that  it  Is  difficult  to  tell  when  a  Stoneley-wave  begins  on 
a  seismogram.  Furthermore,  the  equations  relating  the  Stoneley  wave 
velocity  to  the  shear-wave  velocity  assume  a  liquid  layer  over  a 
homogeneous,  isotropic  solid  medium  which  extends  to  an  indefinitely 
great  distance.  These  conditions  for  the  solid  medium  prove  to  be  an 
ideal  case  that  is  not  seen  In  earth  sediments. 

It  appears  that  one  does  not  have  to  resort  to  determining  the 
shear-wave  velocity  Indirectly  using  Stoneley  waves  or  a  vlscoelastometer. 

By  using  a  specially  designed  SH  source  with  a  correctly  oriented 
horizontal  geophone  and  a  vertical  geophone,  both  the  shear-wave 
velocity  and  compress lonal-wave  velocity  may  be  obtained  from  their 
respective  direct  waves.  (The  horizontal  geophone  is  oriented  to 
produce  a  large  response  to  the  direct  S„  wave  and  the  vertical  geophon. 
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records  the  particle  motion  of  the  direct  P  wave).  This  is  the  method 
investigated  in  this  study. 

The  attenuation  of  the  shear  wave  and  the  compress ional  wave  may 
be  determined  using  the  technique  of  Fast  Fourier  Transforming,  based 
on  a  Cooley-Tukey  algorithm  (Cooley,  1967).  In  this  technique,  one 
obtains  the  amplitude  soectral  density  from  the  analog  voltage  and  then 
squares  this  value.  After  converting  the  resulting  energy  to  decibels, 
energy  versus  frequency  Is  plotted,  and  spectra  at  different  geophones 
are  compared.  For  a  certain  frequency,  the  energy  level  difference 
divided  by  the  geophone  separation  Is  the  ttenuatlon  In  decibels  pet- 
length  unit. 

Density  values  from  laboratory  measurements  on  core  samples  have 
been  obtained  at  the  test  facility  In  San  Francisco  Bay  (section  3.2.8). 
These  values  will  be  used  In  the  equations  to  determine  the  various 
sediment  properties.  Because  the  primary  purpose  of  developing  an  In 
situ  seismic  method  is  to  eliminate  measurements  on  disturbed  core 
samples,  consideration  has  been  given  to  the  measurement  of  the  bulk 
density  in  situ  using  a  nuclear  gauge.  The  use  of  garrma  rays  to  obtain 
bulk  density  in  marine  sediments  has  been  successfully  accomplished  by 
Preiss  (1968)  and  Keller  (1965), 

3.1.3  Viscoelastic  Model 

The  seismic  parameters  to  be  measured  are  the  compress lonal - 
wave  velocity  (Vp) ,  the  shear-wave  velocity  (V5),  the  compress lonal - 
wave  attenuation  coefficient  (ap) ,  the  shear-wave  attenuation  coefficient 
^as),  and  the  bulk  density  (p  ).  These  properties  were  chosen  because 
if  a  viscoelastic  solid  model  is  used  as  a  model  for  the  ocean  floor. 


6k 


then  density,  and  at  least  four  other  physical  properties,  must  be 
known  to  determine  various  other  sediment  properties.  A  viscoelastic 
model  is  favored  over  a  Hookean  model  because  the  Hookean  model  does 
not  explain  the  attenuation  of  seismic  waves  in  marine  sediments 
(Hami I  ton  et  al ,  1970) . 

To  account  for  this  attenuation,  Bucher  0964),  Hamilton,  et  al  (1970), 
and  White  (1965)  replaced  the  rigidity  (p  ),  incompressibility  (K) ,  and 
Lame's  constant  (A)  by  the  complex  Lame  constants  v'Ar  +  ipj), 

( KR+  1  KI )  and  (*n+,*i)  ,n  Lhe  equations  of  elasticity.  However, 
the  linear  viscoelastic  model  for  a  homogeneous  Isotropic  medium 

discussed  by  Flugge  (196?)  and  Fung  <1967)  appears  to  be  the  best  model 
for  the  sea  floor.  The  pertinent  equations  which  result  from  the  theory 
of  Flugge  and  Fung  are  listed  In  section  4.3.2  of  the  appendix. 

Section  4.3.1  summarizes  the  notation  used,  and  section  4.3.3  shows 

the  relationships  between  the  measured  sediment  properties  and  various 
other  proper t I 9S. 


3.1.4  Time-Distance  Curves 

For  the  in  situ  seismic  method  described  In  this  report,  the 
shear-wave  velocity  is  found  from  a  particular  wave  arrival  containing 
information  about  the  shear  wave.  To  determine  the  best  wave  arrival 
from  which  the  shear  wave  velocity  can  be  found,  several  time-distance 
curves  were  constructed  for  a  simple  three-layer  model  where  three 
source-receiver  combinations  were  considered.  The  source-receiver 
combinations  were  an  Impulsive  source  In  the  water  with  the  receiver 
on  the  bottom;  an  Impulsive  source  and  receiver  both  on  the  bottom; 
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and  a  SH-source  and  receiver  both  on  the  bottom  The  method  used  in 
constructing  these  t ime-d I  stance  curves  is  outlined  by  Griffiths  and 
King  (1965).  The  method  consisted  of  calculating  the  time-distance 
relationship  for  each  wave  arrival  and  drawing  its  curve  on  a  time- 
distance  graph.  Multiple  reflections  and  surface  waves  were  neglected 
in  the  calculations.  Reasonable  values  were  chosen  for  the  compress ionai 
wave  velocity  (Vp),  shear-wave  velocity  (Vs),  and  depths  (d)  of  the 
layers.  Figures  3-i ,  3-2  and  3-3  are  examples  of  each  source-receiver 
combination.  In  the  graphs,  the  symbol  "p"  represents  the  compress ionai 
wave,  "Sy"  represents  the  shear  wave  polarized  in  such  a  manner  that  the 
particle  motion  is  in  a  vertical  plane,  and  "S^"  represents  the  shear 
wave  polarized  in  the  horizontal  plane. 

For  an  Impulsive  source  In  the  water  and  receiver  on  the  bottom,  It 
was  concluded  that  wave  arrivals  containing  Information  about  the  shear- 
wave  velocity  would  be  hidden  by  various  reflected  and  refracted  waves 
or  would  have  amplitudes  extremely  small  compared  to  the  expected 
amplitude  of  the  direct  P-wave.  Amplitude  calculations  were  accomplished 
using  the  equations  of  Ewing  (1957). 

For  the  source  and  receiver  on  the  bottom,  it  h*as  concluded  that 
the  direct  Sv-w ave  could  be  theoretically  identified  on  a  seismogram  if 
a  vertical  geophone  is  employed.  All  the  other  wave  arrivals  from 
which  the  shear-wave  velocity  could  be  obtained  were  found  to  arrive 
too  late  and,  consequently,  would  be  hidden  by  various  compress ionai 
waves. 
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Figure  3-1.  Theoretical  time-distance  graph  for  source  on 
water  surface  and  receiver  on  bottom. 
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Figure  3-2.  Theoretical  time-distance  graph  for  both  source 

and  receiver  on  bottom. 
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Figure  3-3.  Theoretical  time-dietanae  graph  for  Su— 
source  and  receiver  both  on  bottom.  H 
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particle  motion  produced  by  the  direct  S„-„ave  is  vertical  and 
perpendicuiar  to  its  line  of  propane, ion.  Thus,  if  a  vertica,  geophone 
is  implanted  in  the  bottom  with  its  sensitive  axis  parallel  to  the 
ground  motion  produced  by  the  direct  S„W.  the  fIrst  ,rrival  on  . 

smogram  should  be  the  direct  Sg-wave.  This  orientation  of  a  verticai 

geophone  to  produce  a  large  response  to  the  direct  Sg-wave  will  be 
referred  to  as  the  “Sy-mode". 

The  compress  iona  I -wave  velocity  for  a  source  and  receiver  on  the 
bottom  may  be  obtained  in  a  similar  manner.  Because  the  particle  motion 
of  a  compression,!  wave  Is  in  the  direction  of  propagation  (I... 
a  horizontal  geophone  with  Its  sensitive  axis  directed  toward  the  sourc,’ 
would  produce  a  large  response  to  the  direct  P-wave.  This  orientation 
of  a  horizontal  geophone  will  be  referred  to  as  the  "P-mode". 

Therefore,  one  of  the  proposed  seismic  methods  was  the  orientation 
of  a  horizontal  geophon.  In  the  P-mod.  and  a  vertical  geophone  In  ,h. 
dg-mode.  The  compress  Iona  I -wave  velocity  and  the  shear-wav.  velocity 
would  be  determined  from  the  direct  P-wave  and  direct  Sg-wave,  respec¬ 
tively.  This  proposed  seismic  method  is  shown  In  figure  3-ba. 

For  the  SH-source  and  receiver  on  the  bottom,  it  was  concluded 

that  the  direct  SH-„ave  was  the  only  wave  that  could  be  efficiently 

identified  on  a  seismogram.  A  horizontal  geophone  could  be  oriented 

*0  that  its  sensitive  axis  is  parallel  to  ,h.  ground  motion  of  the  direct 

s„-wave  (the  ground  motion  being  horizontal  and  perpendicular  to  the 

'ine  of  propagation).  This  orientation  of  a  horizontal  geophone  will 

b.  referred  to  as  the  'Vmode".  Using  the  S^de,  ,h,  first  arrlvai 
should  be  the  direct  S^-wave. 
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Honjontol  Geophone  in  P  -  Mode 
(  Sensitive  Axis  Shown  ) 


a  First  Method 


Figure  o^4. 


Two  proposed  seismic  methods. 
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The  compress  Iona  I -wave  velocity  can  be  obtained  with  a  horizontal 
geophone  in  the  P-mode.  Therefore,  the  second  proposed  method  employed 
the  use  of  a  SH-source  and  two  horizontal  geophones--one  oriented  in  t^.e 
P-mode  and  the  other  oriented  in  the  S^-mode  (figure  3-^b). 

Experiments  were  conducted  to  determine  which  of  the  two  proposed 
seismic  methods  would  be  the  most  feasible  In  marine  sediments. 

3.1.5  Description  of  Experiments 

Experiments  to  determine  the  shear-wave  and  compress ionai -wave 
velocities  were  conducted  In  three  environments.  The  first  environment 
was  a  sand  box  containing  unconsolidated,  moist  sand;  the  second  was  a 
beach  near  the  shoreline  where  the  sediment  was  sand  and  gravel,  and 
the  last  was  In  3  feet  of  water  at  a  beach  on  San  Franc 'sco  Bay  where 
the  sediment  was  also  sand  and  gravel. 

Sand  Box  Tests. 

The  first  experiments,  conducted  in  a  large  sand  box,  were  to 
determine  which  of  the  two  proposed  seismic  methods  would  be  best  to 
obtain  the  shear-wave  velocity.  A  sand  box  was  chosen  because  It  was 
located  in  a  building  where  equipment  and  facilities  were  readily 
accessible.  The  box  was  15  feet  square  and  overlaid  a  concrete  bottom 
at  a  depth  of  3  feet.  Although  the  sand  was  moist  and  very  unconsolidated. 
It  was  unsaturated  so  only  a  rough  idea  of  the  feasibility  of  the  proposed 
seismic  methods  was  determined. 

A  12-pound  sledge  hammer  and  an  impacting  plate  (8  Inches  x  8 
inches)  were  used  as  a  source  of  energy  to  produce,  in  particular, 

Sy-waves  and  P-waves.  The  hammer  contained  an  inertia  switch  so  that 
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the  hammer  impact  would  trigger  the  recording  devices.  The  Su-source 

H 

consisted  of  a  carpenter's  hammer  and  a  metal  stake.  The  Su-source 

n 

and  the  method  of  production  of  S^waves  are  shown  In  figure  3-5.  The 
recorder  was  triggered  as  the  triggering  circuit  was  closed  when  the 
hammer  forced  the  flexible  metal  strip  against  the  base  of  the  stake. 

The  horizontal  geophone  used  in  detecting  the  ground  motion  was  a 
Hall  Sears  geophone  with  a  natural  frequency  of  4  Hz.  The  vertical 
geophone  employed  was  a  Hail  Sears  Model  HS-5-LI  which  has  a  natural 
frequency  of  14  Hz.  These  phones  were  placed  firmly  in  the  sand.  Then 
1-foot  Intervals  were  marked  off  with  a  tap  measure.  From  statlon-to- 
stat ion,  energy  was  transmitted  Into  the  ground  using  the  seismic 
sources,  and  the  resulting  ground  motion  detected  by  the  geophones 
was  recorded  for  the  various  modes  of  operation. 

Beach  Tests. 

Because  satisfactory  results  were  obtained  in  a  moist,  unconsolidated 
sediment  in  identifying  the  arrival  of  the  direct  P-wave  and  direct 
SH-wave,  experiments  were  conducted  next  in  a  slightly  undersaturated 
sediment  employing  two  horizontal  geophones  (one  in  the  Su-mode  and 
the  other  In  the  P-mode).  A  sand  and  gravel  beach  was  chosen  as  the 
field  site.  Using  the  FS-3  refraction  recorder  (sec.  3. 1.6. I),  experi¬ 
mental  runs  were  conducted  on  the  beach  near  the  shoreline.  For  the 

SH-mode  the  SH-source  was  used  while  the  hammer  and  plate  were  used 
for  the  P-mode. 
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Shallow  Water  Tests. 


In  order  to  calculate  the  shear-wave  velocity  and  the  compress ional- 
wave  velocity  In  an  actual  marine  sediment,  a  SH-source  had  to  be  designed 
for  underwater  use.  Several  sources  were  investigated  with  consideration 
given  to  their  energy,  capability  of  producing  S^waves,  and  ease  of 
operation.  Some  of  the  sources  Investigated  were  blasting  caps, 
capacitors,  air  guns,  and  regular  rifle  cartridges.  The  source  finally 
decided  upon  was  a  small  3-cublc-Jnch  air  gun  with  the  capability  of 
being  reduced  to  I -cub ic- inch. 

To  produce  a  large  SH-wave,  a  probe  was  designed  to  be  attached  to 
the  air  gun  (figure  3-6).  As  the  probe  ts  Implanted  Into  sediment, 
water  and  a  little  sediment  will  enter  the  probe.  When  the  air  gun  is 
fired,  the  sudden  expansion  of  the  air  will  force  the  water  and  sediment 
out  the  side  of  the  probe.  The  force  of  this  stream  of  water  and 
sediment  will  cause  a  rotational  movement,  generating  a  large  Sn-wave 
perpendicular  to  the  rotational  movement. 

At  a  beach  on  San  Francisco  Bay  the  air  gun,  with  the  probe  attach¬ 
ment,  was  tested  for  its  effectiveness  of  producing  SM-waves.  The 
field  site  was  In  3  feet  of  water.  The  exit  port  on  the  probe  was 
implanted  by  hand  to  a  depth  of  I  foot  Into  the  sediment.  Two  geophone 
units  (each  contained  a  vertical  geophone  and  horizontal  geophone)  were 
implanted  at  a  distance  of  3  feet  and  5  feet  from  the  air  gun.  The 
horizontal  geophone  was  oriented  in  the  S^-mode.  The  four-trace  oscil¬ 
loscope,  with  camera  attachment,  was  used  for  recording.  The  oscilloscope 
was  triggered  off  the  battery  which  triggered  the  solenoid  of  the  air  gun. 
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Unfortunately,  the  piston,  when  returning  to  its  initial  position 
after  firing  the  air  gun,  sucked  sediment  into  the  probe.  The  probe  of 
the  air  gun  eventually  became  obstructed  with  sediment.  This  limited 
the  number  of  shots  for  every  implant  of  the  probe  (a  maximum  of  three 
shots  achieved).  Sometimes  the  probe  became  obstructed  just  from 
inserting  it  into  the  sediment.  Consequently,  a  little  data  was 
accumulated. 

3.1.6  Data  Acquisition  and  Results 
Refraction  Unit. 

One  method  of  recording  the  wave  arrivals  was  the  use  of  a 
refraction  unit  (Huntec  model  FS-3  Portable  Facsimile  Seismograph), 
a  single-channel  time-distance  plotting  Instrument.  The  record. 

Produced  on  electro-sensitive  paper.  Is  In  the  form  of  short  dashes 
which  represent  the  positive  zero  crossings  of  each  cycle  of  the  wave 
arrivals. 

To  determine  the  feasibility  of  the  proposed  seismic  method 
employing  ’  vertical  geophone  In  the  Sy-mode  and  a  horizontal  geophone 
In  the  P-mode,  experimental  runs  In  the  sand  box  were  conducted  first 
in  the  P-mode  and  then  In  the  Sy-mode  because  the  refraction  unit  Is  a 
single-channel  Instrument.  The  source  employed  was  the  sledge  hammer 
and  impacting  plate.  Figure  3-7a  shows  one  of  the  time-distance  records 
produced  with  the  horizontal  geophone  in  the  P-mode.  The  line  drawn  In 
the  record  represents  the  direct  P-wave.  The  results  for  the  compres- 
sionai-wave  velocity  are  listed  In  table  3-i  for  three  experimental  runs. 
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Figure  3-7 


Time-distance  records  from  unconsolidated  sand 
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Figure  3-7b  is  an  example  of  a  record  produced  with  the  vertical  geophone 
in  the  Sy-mode.  The  velocities  obtained  from  the  inverse  slope  of  the 
first  arrival  for  three  experimental  runs  gave  the  velocity  of  the 
compress  ionai  wave.  The  results  are  listed  in  Table  3-i.  Therefore, 
with  the  refraction  unit,  it  can  be  concluded  that  the  shear-wave 
velocity  cannot  be  obtained  from  the  first  arrivals  with  the  geophone 
oriented  in  the  Sy-mode.  With  difficulty,  It  was  possible  to  draw  the 
curve  that  appeared  to  be  the  direct  Sy-wave  on  some  of  the  records. 


Table  3-1.— Wave  velocities  in  unconsolidated  sand 
using  refraction  unit 


wave 

Compress ionai - 
veloclty(ft/sec) 

Shear-wave 

veiocity(ft/sec) 

First  proposed 

seismic  method 

P-mode: 

6  i  0 

61 0 

63 

- 

Sy-mode: 

665 

610 

6f0 

- 

Second  proposed 

seismic  method 

P-mode: 

- 

- 

S^-mode: 

324 

- 

308 

334 

Next,  the  determination  of  the  shear-wave  velocity  and  compressionai- 
wave  velocity,  using  the  second  proposed  method,  was  conducted  in  the 
sand  box.  With  the  SH-source,  records  were  produced  with  the  horizontal 
geophone  oriented  first  in  the  P-mode  and  then  in  the  SH-mode. 
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Figure  3-8a  shows  an  example  of  the  records  produced  with  the  horizontal 
geophone  oriented  in  the  P-mode.  No  direct  P-waves  were  present  on  the 
records,  revealing  that  the  SH-source  does  not  produce  direct  P-waves  of 
sufficient  energy  to  trigger  the  refraction  unit.  However,  the  direct 
Sn-wave  was  seen  with  the  horizontal  geophone  oriented  in  the  Sn-mode. 
Figure  3-8b  shows  one  of  the  records  and  Table  3-1  ilsts  the  calculated 
vaiues  for  three  experimental  runs.  From  the  beach  tests,  both  the 
direct  S^-wave  and  the  direct  P-wave  were  identifiable  on  the  records. 
Samples  of  the  resuits  are  shown  on  figure  3-9. 

Duai  Beam  Oscilloscope  with  Camera. 

The  second  type  of  recording  instrument  used  was  a  duai  beam 
osciiioscope  (Tektronix  type  555).  Because  this  osciiioscope  has  two 
independent  inputs,  the  outputs  from  two  geophones  were  traced  simul¬ 
taneously  on  the  screen.  To  record  the  wave  arrivals,  a  camera  was 
attached  to  the  osciiioscope  screen.  Polaroid  type  47  photographic 
fiim  (3000  ASA  Speed)  was  used. 

The  same  operational  procedures  were  employed  to  determine  the 
feasibility  of  the  two  proposed  seismic  methods  as  were  conducted  with 
the  refraction  unit.  First,  the  seismic  method  with  a  vertical  geophone 
in  the  Sv-mode  and  a  horizontal  geophone  in  the  P-mode  was  tested  in 
unconsolidated  sand  using  the  siedge  hammer  and  impacting  plate  as  the 
source.  An  exampie  of  a  typicai  picture  obtained  is  shown  in  figure 
3-i0a.  A  time-distance  graph  constructed  from  first  arrivals  from 
severai  pictures  is  shown  in  figure  3-ii  and  reveals  that  the  first 
arrivai  at  both  geophones  is  the  direct  P-wave.  The  results  of  two 
experimental  runs  are  iisted  in  Table  3-2. 
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a.  Output  of  vertical  geophone  in  Sv~mode  (  trace  I  )  and  of 
horizontal  geo-'hone  in  P  -  mode  (  trace  2  )  -  C  sledge 
hammer  and  plate  source,  source  -  receiver  separation  =  4  feet.  ] 


d.  Output  of  vertical  geophone  in  Sv  mode  (trace  I)  and  of 
horizontal  geophone  in  SH -mode  (trace  2)-  [  Sh  source; 
source-receiver  separation  =  I  foot.  ] 

Figure  3-10.  Geophone  records  -  unconsolidated  sand. 
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Figure  3-11.  Time-distance  plots  -  unconsolidated  sand. 

84 


The  second  operation  consisted  of  orienting  the  horizontal  geophone 
in  the  Su-mode.  Because  only  one  horizontal  geophone  was  available, 
experimental  runs  could  not  be  conducted  with  a  horizontal  geophone  in 
the  SH-mode  and  a  horizontal  geophone  in  the  P-mode  (i,e.,  the  second 
proposed  seismic  method).  However,  previous  results  have  shown  that 
the  comp  res sional -wave  velocity  can  be  obtained  from  the  first  arrival 
produced  by  a  vertical  geophone  in  the  Sy-mode;  and  therefore,  a  vertical 
geophone  was  used.  Figure  3- 10b  shows  one  of  the  pictures  taken  using 
the  Su-source.  Plots  of  first  arrivals  on  a  time-distance  graph  for 
two  experimental  runs  are  shown  In  figure  3  —  i 2 .  The  graph  shows  that 
both  the  direct  P-wave  and  the  direct  S^-wave  were  obtalnabie,  and 
their  calculated  velocities  are  listed  in  Table  3-2. 


Table  3-2.— Results  In  unconsolidated  sand 
using  oscilloscope 


Compress lonal- 
veloclty  (ft/sec) 

Shear 

velocity  (ft/sec) 

First  proposed  — 

P-mode: 

580 

-- 

Seismic  method 

Sy-mode’ 

570 

Second  proposal  — 

iy-mode: 

560 

340 

Seismic  method 

SH-mode: 

Additional  experiments  were  conducted  In  unconsolidated  sand  using 
the  dual  beam  oscilloscope.  However,  the  oscilloscope  was  converted  to 
four  traces  by  use  f  a  Tektronix  lA4plug-!n  (I.e.,  the  oscilloscope 
was  converted  to  a  four-channel  recorder).  Furthermore,  four  geophone 
units  were  employed.  Each  unit  consisted  of  a  vertical  geophone  and  a 
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horizontal  geophone  (Geospace  HS-5  Model  K).  These  units  were  implanted 
in  a  line  with  the  horizontal  geophone  oriented  in  the  S^-mode.  A 
typical  picture  for  the  four  horizontal  geophones  oriented  in  the 
S^-mode  is  shown  in  figure  3-13 a.  The  first  arrival  time  of  nach  of 
the  four  traces  was  determined  and  then  plotted  on  a  time-distance 
graph  (figure  3-14).  The  inverse  slope  of  the  line  obtained  from 
these  four  points  gave  the  shear-wave  velocity  from  a  single  picture. 

The  graph  shows  excellent  results  for  five  pictures  taken. 

The  experiment  was  repeated  using  the  four  vertical  geophones  and 
the  hammer  and  plate  source.  A  typical  record  is  shown  in  figure  3-l2b, 
and  the  results  of  four  pictures  are  plotted  on  a  time-distance  graph 
in  figure  3-15.  The  compress ional  velocities  obtained  are  listed  In 
the  graph. 

For  shallow  water  tests,  one  of  the  problems  encountered  In  the 
oscilloscope  recording  was  reading  data  values  that  were  burled  due  to 
too  large  a  time  scale  on  the  oscilloscope.  It  was  difficult  to  measure 
the  arrival  time  difference  of  the  direct  P-wave  between  the  two  vertical 
geophcnes.  Figure  3-16  shows  two  pictures  that  were  taken.  By  measuring 
the  difference  In  arrival  time  of  the  direct  P-wave  and  direct  S^-wave 
between  their  respective  geopnones  and  dividing  by  the  geophone 
separation  (2  feet),  a  compress ional -wav.'  velocity  of  5,000  ft/sec 
was  calculated  from  both  pictures,  and  a  shear-wave  velocity  of  340 
ft/sec  and  242  ft/sec  was  calculated  from  picture  1  and  2  of  figure 
3-1 6  respectively. 
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a.  Wave  arrivals  from  four  horizontal  geophones  in  Sh  *  mode 
(source-receiver  separation  distance  listed  cn  traces.  ) 


b.  Wave  arrivals  from  four  vertical  geophones  in  Sv  "  mode 
(  source -receiver  separation  distance  listed  on  traces.  ) 


Figure  3-13.  Geophone  records  -  unconsolidated  sand. 
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Figure  3-14.  Time-distance  plots  -  unconsolidated  sand. 
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3-1.7  Conclusions  and  Recommendations 

Of  the  two  proposed  seismic  methods  studied,  only  one  proved 
to  be  feasible  in  determining  the  compress ional -wave  velocity  and  shear- 
wave  velocity.  The  method  consisting  of  a  vertical  geophone  oriented  to 
produce  a  large  response  to  the  direct  Sy-wave  and  of  a  horizontal 
geophone  to  produce  a  large  response  to  the  direct  P-wave  was  shown  to 
give  only  the  compress ional -wave  velocity  and  not  the  shear-wave 
velocity  from  first  arrivals.  This  was  due  to  the  fact  that  either 
the  vertical  geophone  was  not  implanted  perfectly  normal  to  the  surface 
(if  one  assumes  that  the  direct  P-wave  travels  on  the  surface),  or 
velocity  Increases  with  depth  in  the  sediment,  which  would  result 
in  the  direct  P-wave  making  an  arcuate  path  to  the  surface.  Both 
would  cause  the  ground  motion  of  the  direct  P-wave  to  be  detected  by 
the  vertical  geophone. 

The  second  proposed  seismic  method  was  shown  to  give  both  the 
shear-wave  velocity  and  compress iona I -wave  velocity  from  first  arrivals 
in  unconsolidated  sand,  on  a  beach  of  slightly  undersaturated  sand  and 
gravel,  and  in  San  Francisco  Bay  in  saturated  sand  and  gravel.  Using 
an  SH-source,  the  direct  SH-wave  was  seen  in  the  records  with  a 
horizontal  geophone  oriented  to  produce  a  large  response  to  the  direct 
SH-wave.  For  this  method,  the  direct  P-wave  can  be  seen  either  with  a 
horizontal  geophone  oriented  to  produce  a  large  response  to  the  direct 
P-wave  or  with  a  vertical  geopho’ e. 

As  may  be  seen  from  the  experiments  conducted,  the  primary  goal  to 
date  has  been  the  determination  of  a  seismic  method  to  obtain  the  shear- 
wave  velocity.  It  Is  believed  that  by  using  the  air  gun  with  the  probe 
attached,  SH-waves  will  be  generated  with  sufficient  energy  to  be  detected. 
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After  satisfactory  results  have  been  obtained  at  the  beach  using 
the  air  gun,  an  acoustically  dead  framework,  on  which  will  be  attached 
the  four  geophone  units  and  the  air  gun,  will  be  designed  and  constructed 
A  sketch  of  the  framework  is  shown  in  figure  3-17.  Upon  completion  of 
the  framework,  the  stations  in  San  Francisco  Bay,  where  core  samples 
were  taken  (section  3.2.8),  will  be  reoccupied.  After  determining 
the  compress ional-wave  velocity,  shear-wave  velocity,  and  their 
attenuation  coefficients,  and  using  the  value  for  bulk  density  deter¬ 
mined  from  the  core  samples,  the  sediment  properties  in  Table  3-3  will 
be  calculated.  Furthermore,  an  attempt  will  be  made  to  find  empirical 
relationships  from  which  porosity,  shear  strength,  bulk  density,  and 
mean  grain  size  can  be  predicted  from  the  In  situ  measurements  of  the 
compress Ional-wave  velocity  and/or  the  shear-wave  velocity.  These 
four  sediment  properties  have  been  related  to  the  compress ional-wave 
velocity  for  various  environments  by  Hamilton  (1970b),  Horn  (1968), 
Schreiber  (1968),  and  Lewis  (1971). 

Because  the  undisturbed  bulk  density  Is  preferable  for  calculating 
the  sediment  properties  for  a  viscoelastic  solid  model,  consideration 
has  been  given  to  the  construction  of  a  nuciear  density  probe  to  be 
attached  to  the  framework.  Cesium  137  with  a  strength  of  5  to  10 
mil  1 icur ies  seems  to  be  suitable  radioisotope,  it  has  a  half  life  of 
approximately  30  years  and  emits  a  gamma  photon  with  a  discrete  energy 
of  .662  mev.  A  scintillation  counter  would  be  used  to  detect  the 
electrons  emitted  by  the  sediment.  The  separation  distance  would  be 
approximately  I  foot. 
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3.2  Classification  of  Sea-Floor  Sediments  by  a  Ship  Underway 
3.2.1  Concept 

Fhe  geologic  nature  of  the  sea  floor  hcs  long  been  known  to  have 
certain  intrinsic  physical  properties  amenable  to  measurement  by  remote 
sensors  as  well  as  in  situ  devices.  Acoustic  soundings  of  the  bottom 
and  subbottom  with  seismic  energy  in  the  form  of  short  pulses  and  re¬ 
cording  the  echo  or  reflected  energy  from  acoustically  reactive  horizons 
has  become  a  popular  method  employed  by  geophysicists  to  map  the  vari¬ 
ability  of  the  earth's  geologic  structure.  Although  it  has  been  theo¬ 
retically  recognized  that  the  method  has  application  for  determining 
mass  physical  properties  of  sediments  and  rock,  the  utilization  has  not 
developed.  Requirements  for  this  kind  of  data  were  minimal  until  ocean 
exploration  began.  The  need  for  tea  floor  engineering  properties  and 
better  knowledge  of  sediment  distribution  has  prompted  new  interest. 

Two  obstacles  still  remain  and  have  limited  widespread  application  of 
the  methods.  The  first  concerns  the  time  consuming  detail  required  in 
data  processing  and  the  second  lies  in  the  uncertainty  involved  in 
relating  the  attenuation  or  loss  of  seismic  energy  noted  in  continuous 
seismic  profiling  to  parameters  depicting  the  engineering  and  physical 
properties  of  the  materials.  The  advent  of  the  digital  computer  has 
overcome  the  first  obstacle,  but  many  more  field  exercises  must  be 
conducted  to  take  measurements  before  a  thorough  understanding  of  the 
uncertainties  of  the  second  obstacle  is  possible.  The  work  described 
herein  entails  one  such  exercise. 
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5-2.2  Procedure 

A  si.pl,,  straightforward  approach  to  making  acoustic  measurements 

°f  ^  “  f,°°r  ^  d'SCribed-  i*  involves  ,h,  continuous 

.app.ng  of  reflection  coefficients  by  coppering  the  peak  amplitude  of 

*"C  tranSmi"ed  PUl!'  «*“*“>  “  «-  P-k  amplitude  of  the  echo  or 
tef.ec.ed  puls,  (response,.  this  investigation  the  acoustic-ref lee 
•tvi.y  measurement  is  made  at  norma,  incidence  and  with  both  a  wide-band 

°'S  WZ  <W‘th  Ce",er  tTf^y  at  2.5  KHz)  sparker  source  and  a  discrete 
A  semi -automated  measurement  system  was  designed  and 
constructed  to  facilitate  the  acquisitio„  of  acoustic  data  on  shipboard 
Continuous  recording  of  seismic  events  depicting  the  i^ulse-response 
functions  for  each  of  the  two  sources  and  their  respective  receivers  is 
obtained.  The  peak  amplitude  measurement,  are  taken  from  the  zen  Un, 
of  the  pressure  wave  train,  measuring  only  posltive  exc(lrsion  ^ 
line. 

data  presented  in  this  report  was  developed  without  the  aid  of 
digital  computer  processing.  Analog  measurements  were  reduced  manually 
from  photographic  reproduction  of  data  acquired  on  analog  magnetic  tape, 

»  was  felt  that  the  initial  investigation  shold  pursue  this  course  to 
gain  better  familiarization  with  the  sign.l/nois,  variations  and  under- 

standing  the  physical  (geometricall  nrohu., 

vg  o metrical;  problems  encountered  in  making 

acoustic  measurements  of  this  type.  Future  plans  will  include  machine 

processing  which  will  undoubtedly  enhance  the  precision  and  handling  of 
the  voluminous  quantity  of  data. 
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3.2.3  Previous  Investigations 

The  first  analysis  of  acoustic  transmission  in  seawater  began 
aoring  horld  bar  II  when  several  investigators  looked  into  submarine 
warfare  applications  11946).  Slant  angle  measurements  at  various  f re¬ 
gencies  and  distances  provided  the  first  data  on  the  sea  floor  geological 
environment,  attention  was  focussed  on  two  aspects  relating  to  wave 
propagation  phenomena  and  sediments  per  se.  In  the  period  following  the 
war,  investigations  by  Liebermann  (1948).  Urick  (1984),  Uriel  and  Sal.ng 
(11162),  Mackenzie  (1960),  McKinney  and  Anderson  (1964,  and  Jones,  et  al 
0964)  were  directed  to  the  characteristics  of  sound  transmission  and 
reflection.  Attention  was  not  given  to  sea-floor  sediments  and  their 
mass  physical  properties  until  the  mid  1950s  when  investigations  by 
Hamilton  (1956),  Shumway  (1960),  Sarmiento  and  Kirby  (1962),  Richards 
0962,  and  Nafe  and  Drake  (1963,  commenced.  Actual  field  measurements 
of  reflectivity  and  correlations  to  sediment  type,  by  a  ship  undetway, 
were  made  by  Taylor-Smith  (1968),  Taylor-Smith  and  Li  (1966)  and 

Bres.au  (1964  and  1967)  in  the  mid  1960s. 

■me  material  presented  in  this  section  is  another  effort  to  obtain 
fteld  measurements,  but  with  a  practical  application  in  mind.  The  work 
generally  follows  the  recommendations  of  Breslau  (196?  and  perso 
communications)  that  more  investigations  and  measurements  are  needed  over 
all  water  covered  areas.  This  effort  for  puiposes  of  satisfying  both 
the  objectives  of  the  contractor  (ARPA)  and  NOAA  is  limited  to  the 
Continental  Shelf  water  depths.  It  is  felt  the  practical  application 
for  sea-floor  tunneling  and  excavation  for  foundations  engineering, 
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mining  (dredging)  and  defense  purposes  necessitate  concentrating 
attention  on  this  zone.  Several  technologic  advances  since  the  work 
performed  by  Breslau  (1964  and  1967)  are  incorporated  in  the  electronic 
systems  id  in  this  investigation. 

3.2.4  Theory  Developed 

A  simplified  approach  to  a  very  compU'  physical  process  involving 
wave  propagation,  attenuation  and  absorption  has  been  used  in  develop¬ 
ing  the  model  of  the  sea  floor-water  interface.  The  sea  floor  is 
considered  as  a  plane  interface  between  two  fluids.  Energy  is  focussed 
on  a  zone  of  about  one  square  meter  surface  area.  The  backs cattering 
strength  associated  with  this  surface  area  depends  much  on  the  roughness 
of  the  sea  bottom.  A  smooth  bottom  acts  as  a  mirror  such  that  the  re¬ 
flected  energy  is  concentrated  on  the  specular  direction.  A  rough 
bottom  will  reflect  the  energy  more  evenly  in  all  directions.  It  is 
important  to  note  that  the  terms  "smooth'*  and  "rough"  must  be  considered 
in  relation  to  the  acoustical  wave  length.  Where  the  wave  length  is 
only  a  few  centimeters  it  is  seen  that  a  100%  smooth  bottom  is  rare. 

Even  small  stones,  holes  or  ripple?  will  make  the  bottom  rough  in  the 
sense  considered  here.  When  f.  sound  wave  strikes  the  sea  bottom,  part 
of  the  incident  energy  will  be  reflected.  The  reflected  energy  is 
spread  in  all  directions  with  a  typical  distribution  as  shown  in  figure 
3-18.  The  distribution  is  totally  dependent  on  the  roughness  of  the 
sea  bottom.  Figure  3-19  illustrates  how  the  roughness  of  the  topography 
can  influence  the  bottom  backscattering  strength  as  a  function  of  the 
grazing  angle.  The  angle  of  incidence  is  seen  to  also  play  an  important 
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Figure  3-18.  Seafloor  scattering  -  specular  energy. 
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BOTTOM  BACKSCATTERING 


Figure  3-19.  Typical  curves  for  bottom  baokeoattering 
strength  (After  Bodholt  1969). 
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role  in  determining  the  strength  of  the  echo  signal.  The  experiments  in 
the  following  description  were  all  conducted  with  the  source  and  receiver 
set  up  to  record  normal  incidence.  It  can  be  shown  that  theoretically 
the  most  accurate  measurements  of  bottom  backscattering  strength  should 
be  obtained  with  a  narrow  beam  transducer.  If  the  beam  width  is  5-10° 
between  the  3dB  points  the  bottom  backscattering  strength  will  vary  only 
a  few  dB  or  less  within  the  beam,  of  course  depending  on  the  type  of 
bottom.  The  magnitude  of  the  bottom  sampling  area  will  depend  on  the 
pulse  length  of  the  transmitted  sound.  The  shorter  the  pulse  the  smaller 
the  sampling  area.  Figure  3-20  illustrates  the  beam  width,  pulse  length 
and  area  concept. 

In  the  exercise  conducted  for  this  investigation,  two  sources  are 
employed.  One  is  an  uncontrolled  beam  width  emitting  a  .5  millisecond 
pulse  length  and  the  other  is  a  narrow  beam  (7.5°)  echo  sounder  having 
a  pulse  length  of  .25  millisecond. 

Derivations. 

A  detail  treatment  of  the  mathematical  description  for  acoustic  wave 
propa' ition,  reflection  and  absorption  (bottom  loss)  is  contained  in  the 
appendix  (section  4.4)  to  this  report.  Several  references  were  consulted 
and  the  best  mathematical  expressions  relating  to  the  experiments  con¬ 
ducted  in  this  investigation  were  taken  from  derivations  by  Breslau 
(1968). 

From  the  derivations  the  most  useful  equations  applicable  to  data 
reduction  and  processing  are  represented  by  defining  the  relationship 
between  the  fractional  loss  of  intensity  at  the  sea  floor,  K,  and  the 
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Direction 


Figure  3-20.  Transmitted  sound  pulse 


(Eq.  13) 


Rayleigh  Reflection  Coefficient,  R. 

1 .e.  k  =  R 

where  K  is  defined  on  a  peak  pressure  basis,  thusly 

K  ’  (2D>2  <*■  8) 

where 

PR  =  pressure  of  the  echo  pulse 

Pg  =  pressure  of  the  source  pulse 

D  =  transmission  loss  due  to  spherical  spreading 

e"a2D  =  transmission  loss  due  to  dissipative  attenuation  of 
sound  in  seawater 

Expressing  Eq.  15  as  a  logarithmic  function  to  get  decibel  form  gives 
an  expression  of  the  fundamental  equation  for  bottom  loss,  BL 

i.e.  BL  =  SpL  -  EpL  -  TLS  -  TL^  (Eq.  15) 

where 

BL  =  bottom  loss  (in  dB)  and  is  equal  to  -  10  log  K 

2 

SPl  =  pressure  level  of  the  source  pulse  (in  dB//l  dyne/cm  ) 
and  is  equal  to  20  log  Pg 

EpL  *  pressure  level  of  the  echo  pulse  (in  dB//l  dyne/cm2) 
and  is  equal  to  20  log  PR 

TLg  =  spreading  loss  (in  dB)  and  is  equal  to  20  log  2D 
TL^  *  attenuation  loss  (in  dB)  and  is  equal  to  10  2D  log  e 
A  final  expression  relates  to  bottom  loss,  BL,  to  reflectivity,  R. 

BL  =  20  log  R  (Eq.  17) 

Other  manifestations  of  these  basic  expressions  that  include  measure¬ 
ments  dealing  with  the  energy  versus  time  waveform  of  the  echo  are 
described  in  Breslau's  paper  (1968).  Since  this  investigation  was 
concerned  only  with  peak  pressure  measurements,  the  derivations  relating 
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to  energy  in  the  echo  strength  are  not  included  in  the  appendix. 

Relation  of  Bottom  Loss  to  Mass  Physical  Properties. 

The  mass  characteristics  of  the  bottom  sediment  are  related  to 
bottom  loss  through  the  Rayleigh  Reflection  Coefficient  and  the  acoustic 
impedance  contrast  (Eq.  12)  at  the  interface,  and  it  has  been  shown  by 
Eq.  17  how  the  Rayleigh  Reflection  Coefficient  is  related  to  bottom  loss. 
It  is  possible  to  establish  a  relationship  between  porosity  and  bottom 
loss  through  the  expression  for  impedance  (  Z  )  which  is  dependent  upon 
density  and  compressional  velocity,  both  factors  that  relate  to  the 
porosity  of  natural  sediments. 

In  summary,  several  investigations  have  ascertained  a  linear  re¬ 
lationship  exists  between  porosity  and  density,  in  particular  that 
bottom  loss  as  deduced  from  many  measurements  made  on  sediments  bears  a 
definite  relationship  to  the  porosity  of  the  sediment.  And  furthermore, 
since  the  porosities  of  natural  sediments  are  somewhat  related  to  sedi¬ 
ment  grain  size,  it  is  also  possible  to  establish  a  general  relationship 
between  bottom  loss  and  the  geological  properties  of  sediments. 

For  further  discussion  of  the  sediment  bulk  properties  and  vari¬ 
ations  within  the  physical  makeup  of  the  grains,  interstitial  properties, 
etc.  and  more  specifics  on  how  these  factors  influence  porosity,  the 
reader  is  directed  to  the  papers  by  Nafe  and  Drake  (1963),  Birch  et  al 
(1942),  Hamilton  et  al  (1956),  Shumway  (1960),  Sutton  et  al  (1957)  and 
Breslau  (1964).  The  information  available  on  this  subject  is  too  vo¬ 
luminous  to  be  discussed  here. 
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3.2.5  Description  of  Geophysical  Instrumentation  and  Equipment 

Ihe  instrumentation  used  for  acquisition  and  processing  the  data 
is  discussed  independently  according  to  its  location  and  function.  The 
recording  and  sound  source  equipment  is  all  shipborne,  the  reproducing 
and  processing  equipment  is  located  in  a  shore-based  central  processing 
facility.  The  block  diagrams  given  in  figures  3-21  and  3-22  illustrate 
the  data  flow.  Elements  in  the  data  flow  are  indexed  by  circled  letters 
for  purposes  of  discussing  the  components  and  their  respective  signal 
processing  function. 

Acoustic  Sources  (Shipborne). 

The  acoustic-reflectivity  measurements  incorporated  two  acoustic 
impulse  sources,  an  electrical  type  discharge  generated  by  a  multiple- 
electrode  sparker  sound  source  (figure  3-23)  and  a  standard  fathometer 
transducer.  The  directional  multiple-electrode  sound  source  (DMSS), 
developed  by  MMTC,  is  powered  by  a  capacitor  bank  capable  of  delivering 
a  maximum  output  of  1000  joules.  For  this  experiment,  500  joules  was 
the  operating  power,  measured  at  the  capacitor.  Effective  power  de¬ 
livered  at  the  source  varies;  however,  tests  indicate  at  least  a  60% 
power  loss  occurs  due  to  transmission  loss.  This  is  a  factor  normally 
overlooked  in  rating  power  delivered  in  spark  discharge  seismic  systems. 

A  more  meaningful  value  is  the  acoustic  output  determined  by  calibration 
using  a  standard.  The  measured  average  acoustic  power  output  of  the 
source  is  101.5  dB/V/ubar  referenced  to  1  yard. 

The  DMSS  source  is  a  prototype  test  tool  that  produces  a  broad 
band  frequency  spectrum  which  ranges  from  0  to  5  KHz.  Center  frequency 


?05 


Tektronix 
Type  3150 
Oscilloscope 


1 - - 

I _ 


‘'□IG*  Audio 


Audio 


•'night  Kit 
Duol  Input 
Audio  Amplifier 


j  Wovet< 
r — j  vcg 


Wovetek  Model  III 
VCG  Function 
Generotor 


OTik  tram 

nr 

ClCill^icooi 


n 


Honeywell 
Model  Moo 
Seven-  T;«e  , 
Mognetlc  Tap* 
j  Recorder 
I  Direct 
Reproduce 
Ampllfiere 


Eet  Tcfre  com  rat 


'1NOT£:  Sfgnol  Output 

1  Time  Code 

2  Fothometer  Echoes 

3  Impulse  Hydrophone 

2  Fothometer  Tronsmltter 

3  Sfreomer  Hydrophone 


MMTC 
Vorloble 
Time  Drive 
Unit 


Honeywell  Model  1500a" 
Visicorder  with  5  Type 
_M  -  500Q  Golvonomet(»r« 


| _ Test  Leod 


mucing  equipnent. 


107 


is  about  2.5  KHz.  The  short  pulse  length  and  wide  beam  enable  a  specular 
energy  distribution  yet  retain  a  relatively  high  resolution  for  shallow 
(water  operation)  penetration.  Directionality  of  this  source  has  been 
tested  giving  about  3  to  5  dB  improvement  in  the  vertical  (downward 
projection)  plane  when  the  cylinders  are  in  the  downward  orientation. 

The  multiple  electrode  arrangement  produces  what  is  essentially  a  line 
source.  Bubble  collapse  which  is  known  to  generate  undesirable  secondary 
impulses  is  minimal.  This  was  the  primary  reason  the  DMSS  unit  was  se¬ 
lected  over  single  and  triple  electrode  sources. 

In  conjunction  with  the  DMSS  unit,  a  Raytheon  Model  DE-731  Pre¬ 
cision  Fathometer  and  Depth  Recorder  (figure  3-24)  was  used.  The  41.5 
KHz  discrete  frequency  source  is  a  hull-mounted  barium  titanate  piezo¬ 
electric  type  transducer.  It  has  a  narrow  beam  angle  of  16°  at  -3dB  on 
the  frontal  plane  and  a  43°  beam  in  the  athwart  ships  plane.  The  trans¬ 
mitted  pulse  is  rated  at  600  watts  into  a  150  ohm  load.  In  terms  of 
acoustic  power  delivered  this  is  113  dB/V/ubar  at  1  yard. 

Impulse- Response  Receivers  (Shipbome). 

The  sparker  impulse  (transmitted  peak  pressure)  hydrophone  receiver 
array  (figure  3-25)  consists  of  nine  hydrophones  connected  in  parallel 
whose  output  is  fed  directly  to  a  shipboard  single  channel  amplifier. 

The  array  is  mounted  piggyback  to  the  acoustic  source  unit  during  the 
operations,  thus  representing  as  near  a  direct  coupling  to  the  source  as 
possible  (figure  3-26).  Throughout  the  survey,  the  amplifier  gain  of 
the  piggyback  hydrophone  array  is  set  at  near  unity  to  eliminate  over¬ 
loading  and  signal  distortion  that  could  affect  the  amplitude  measure¬ 
ment. 
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The  sparker  source  with  the  piggybank  hydrophone  is  towed  alongside 
a  main  hydrophone  streamer  at  a  distance  of  about  75  feet  astern  of  the 
ship.  The  main  hydrophone  streamer  receives  the  reflected  response  from 
the  bottom  and  subbottom  seismic  horizons.  The  maximum  side-by-side 
distance  employed  throughout  the  survey  was  about  10  feet,  thereby 
affording  the  reception  of  essentially  a  normally  incident  signal. 

The  average  effective  dynamic  range  (in  field  operation)  of  the 
impulse  hydrophone  is  80  dB.  The  low  end  sensitivity  is  limited  by  random 
electronic  noise.  The  basic  receiving  sensitivity  is  -92.4  dB/V/ubar, 
and  the  maximum  amplifier  gain  is  42  dB. 

The  signals  from  the  sparker  impulse  hydrophone  and  the  main  hydro¬ 
phone  streamer  are  fed  independently  through  the  shipboard  amplifier  and 
then  onto  separate  tracks  of  a  7  channel  analog  magnetic  tape  recorder 
(figure  3-27).  Shipboard  amplification  of  the  impulse  receiver  is  kept, 
near  unity.  The  main  hydrophone  streamer  amplification  varies  according 
to  the  gain  needed  to  boost  the  signal  above  the  noise  threshold.  Con¬ 
ditions  that  can  cause  changes  in  the  gain  are;  1)  large  variation  in 
water  depth,  2)  externally  generated  noise,  or  3)  marginal  reflectivity 
due  to  high  energy  loss  in  soft,  unconsolidated  sediments.  Annotation 
of  the  variations  in  gain  settings  to  the  amplifiers  is  recorded  as 
changes  are  made.  The  data  is  then  normalized  in  processing. 

In  some  instances  the  signal  generated  by  the  number  1  main  hydro¬ 
phone  streamer  is  inverted.  This  is  due  to  processing  the  signal 
through  a  100  Hz  high-pass  filter.  The  graphical  reproduction  of  the 
analog  of  the  wave  train  on  the  visicorder  display  will  therefore  re¬ 
flect  this  inversion.  Whether  it  is  inverted  or  not  is  dependent  upon 
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Figure  3-27.  7-channel 
analog  tape  recorder 
(Ampex  model  SP-300) 
and  receiver-amplifier . 


2  main  hydrophone  array 
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filtering.  The  number  2  main  hydrophone  streamer  signal  is  always  in¬ 
verted. 

Two  main  hydrophone  streamers  were  alternated  throughout  the  survey. 
One  streamer,  designated  as  No.  1,  (figure  3-25)  consisted  of  20  hydro¬ 
phone  elements  divided  into  8  electrical  sections  or  groups  of  either  2 
or  3  elements.  Within  each  group,  the  phones  are  wired  in  parallel.  The 
groups  are  wired  in  series.  The  signals  from  the  elements  are  amplified 
by  preamplifier  modules  in  the  streamer.  The  number  of  groups  activated 
vary  according  to  the  operator's  selection.  Variation  is  controlled  at 
the  shipboard  amplifier  unit.  The  signals  from  all  8  groups  or  fewer 
can  be  composited  and  amplified  at  the  shipboard  control  unit.  The 
average  effective  dynamic  range  of  the  number  1  hydrophone  is  29.2  dB. 

The  low  end  sensitivity  is  limited  by  the  towing  noise  observed  on  the 
records.  The  nominal  operating  dynamic  range  is  46  dB  where  the  low  end 
sensitivity  is  limited  only  by  ambient  noise  generated  by  the  electronics 
and  vibrations.  The  basic  receiving  sensitivity  is  -75.5  dB/V/ubar.  The 
maximum  amplifier  gain  is  52.3  dB  with  the  100  Hz  filter  on. 

A  depth  transducer  is  contained  in  the  number  1  streamer  and  infor¬ 
mation  regarding  steamer  towing  depth  is  read  out  on  a  calibrated  meter 
on  the  front  panel  of  the  shipboard  amplifier  unit. 

The  number  2  main  hydrophone  streamer  (figure  3-28)  consists  of  18 
hydrophone  elements  divided  into  3  electrical  groups  of  6  elements,  each 
wired  in  parallel,  with  the  groups  wired  in  series.  Preamplification  of 
the  signal  and  shipboard  amplification  and  capability  for  group  vari¬ 
ation  is  the  same  as  used  on  the  number  1  streamer.  The  average  effec¬ 
tive  dynamic  range  (in  field  operation)  of  this  hydrophone  is  32  dB. 
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The  low  end  sensitivity  is  limited  in  this  case  by  towing  noise  observed 
on  the  records.  The  nominal  dynamic  range  is  47  dB.  The  low  end  sensi¬ 
tivity  in  this  case  is  limited  only  b}  ambient  electronic  and  vibration 
noise.  The  basic  receiving  sensitivity  is  -76  dB/V/ubar  and  the  maximum 

amplifier  gain  is  46  dB.  The  number  2  streamer  does  not  have  a  depth 
transducer. 

The  receiver  hydrophone  for  the  41.5  KHz  fathometer  transducer  has 
a  dynamic  operating  range  of  60  dB.  This  value  was  not  measured,  but 
was  obtained  from  the  Raytheon  Model  DE-731  Fathometer  Depth  Recorder 
Operation  and  Maintenance  Manual. 

Signal  Processing  (Shipbome). 

lhe  transmitted  impulses  and  return  echoes  are  monitored  at  the 
transducer  input.  The  signals  are  next  fed  to  attenuator  circuits  to 
reduce  the  transmitted  impulses  to  an  acceptable  recording  level.  The 
transmitted  signal  was  attenuated  by  43.2  dB  and  then  fed  directly  to 
the  magnetic  tape  recorder.  The  return  echoes  are  attenuated  by  only 
2.7  dB,  while  the  accompanying  transmitted  signal  was  clipped  by  diodes. 

The  low  amplitude  echoes  are  then  amplified  and  these  signals  fed  to  the 
magnetic  tape  recorder.  The  transmitted  fathometer  impulse  amplifier  is 
an  R-C  coupled  wide  band,  solid  state  amplifier  with  a  40  dB  maximum  gain 
and  a  frequency  response  of  20  Hz  to  500  KHz.  The  amplifier  features 
low  harmonic  distortion,  flat  frequency  response,  and  high  stability  due 
to  the  fact  that  its  variable  gain  is  determined  by  a  ratio  resistor  that 
vanes  the  amount  of  negative  feedback  rather  than  a  nonlinear  transistor. 
All  specifications  for  the  Wide  Band  Amplifier  (Model  1202-3)  were  taken 
from  the  Technical  Manual  of  the  manufacturer,  California  Electronic 


mmmm  »* 


115 


Manufacturing  Co. 

The  Ampex  Model  SP-300  rack  mounted  magnetic  tape  recorder  is  a 
seven  track  recorder  using  1/2  inch,  1.0  mill  tape  on  type  NAB  10-1/2 
inch  reels.  It  is  operated  at  a  tape  speed  of  15  ips,  accurate  to  a 
rated  +_  0.2 %  and  with  flutter  specification  of  less  than  0.15%.  The 
head  assemblies  are  I RIG  compatible.  The  seven  individual  plug-in 
amplifiers  enable  direct  or  FM  record  or  reproduce  capability.  The 
direct  record  mode  of  operation  is  used  in  this  application. 

The  input  signal  range  required  for  a  normal  level  of  recording  is 
specified  at  0.5  to  10.0  volts  rms  with  total  harmonic  distortion  rating 
of  less  than  0.5%  tested  at  1  KHz  and  15  ips.  The  input  impedance  is 
greater  than  10  K  ohms.  The  record  amplifier  can  produce  an  output  to 
the  head  at  least  8  dB  above  the  norma?  level  before  saturation  of  the 
amplifier  occurs.  The  rated  signal  to  noise  ratio  is  30  dB.  At  15  ips, 
the  rated  frequency  response  is  within  +_  3  dB  from  50  Hz  to  40  KHz  & 

During  field  operations,  selected  calibration  signal  levels  are 
recorded  on  each  tape.  At  the  shore  based  playback  center,  these  pre¬ 
recorded  calibration  signals  are  used  to  determine  any  signal  level 
errors  of  both  tape  recorders  and  the  visicorder.  The  errors  are  then 
added  to  the  computed  signal  levels. 

The  inputs  to  the  seven  channels  are  as  follows  (figure  3-21): 

Channel  1:  Voice 

Specifications  quoted  were  obtained  from  SP-300,  FM/Direct 
Recorder/Reproducer  Operations  and  Maintenance  Manual,  Ampex  Part.  No. 
48814-30. 
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Channel  2: 


IVWV  audio 


Channel  3:  Fathometer  transmitted  impulse  signals  (figure  3-29e) 

Channel  4:  Fathometer  echos  (reflected  signals)  (figure  3-29d) 

Channel  5:  Sparker  impulse  signals  from  the  piggyback  hydro¬ 
phone  (figure  3- 29b) 

Channel  6.  Time  code,  IRIG  standard  format  B  (figure  3-29c) 
Channel  7:  Main  streamer  hydrophone  (reflected  signals)  (figure  3 
The  time  reference  is  obtained  from  the  Systron-Donner  Model  8120 
Time  Code  Generator.  The  unit  is  used  to  generate  an  Inter  Range  In¬ 
strumentation  Group  R  (IRIG  B)  time  code  format  in  terms  of  Binary  Coded 
Decimal  (BCD)  hours,  minutes,  and  seconds.  As  a  Time  Code  Generator,  it 
also  has  a  precise  digital  slock,  which  produces  a  serial  code  for  re¬ 
cording  continuous  time  on  the  magnetic  tapes.  On  the  front  panel  of 

the  instrument,  updated  time  is  displayed  by  an  in-line  nixie-tube  read¬ 
out  unit. 

The  time  base  is  a  crystal  controlled  oscillator  with  an  overall 
stability  of  1  part  in  105  per  day.  The  modulated  code  is  generated  on 
a  1  KHz  carrier  with  an  adjustable  amplitude  from  0  to  10  volts  (peak  to 
peak)  and  an  adjustable  modulation  ratio. V  This  output  was  connected 
directly  to  the  Ampex  tape  recorder. 

Signal  Processing  (Shore  Based). 

The  recorded  tapes  are  processed  at  a  shore  based  playback  center 
(figure  3-30).  The  tapes  are  first  played  back  at  15  ips  (recording 

Information^from-Systron-Donner  Time  Code  Generator  Model  8120 
Manual . 
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Figure  3-29.  Signal  ware  forms. 
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Signal  wave  forms 


Figure  3-30.  Shore-base  playback  console. 


speed)  to  listen  to  the  audio  channel  and  note  the  gain  setting  changes 
or  pertinent  information  which  might  affect  the  computation  of  the  signal 
levels. 

The  Honeywell  Model  5600  rack  mounted  magnetic  tape  system  is  a  7 
track  machine  capable  of  reproducing  data  at  seven  selectable,  servo- 
controlled  tape  speeds  from  15/16  ips  through  60  ips.  The  machine  has 
TRIG  compatible  1/2  inch  tape  heads  with  an  optional  voice  channel  on  an 
edge  track.  Flutter,  skew,  and  time  displacement  eiror  are  minimized  by 
the  servo-controlled  tri-capstan  drive  and  independent  reel  servos.  At 
15  ips,  flutter  is  rated  at  0.4%  P-P  and  at  1-7/8  ips  0.7%  P-P.  Tape 
speed  accuracy  is  rated  at  0.15%  of  the  selected  speed. 

The  machine  is  fitted  with  7  direct  reproduce  plug-in  amplifiers 
and  one  direct  plug-in  amplifier.  The  record/reproduce  amplifier's 
rated  frequency  response  at  15  ips  is  100  Hz  to  75  KHz  with  a  30  dB  RMS 
signal  to  RMS  noise  ratio;  at  1-7/8  ips,  50  Hz  to  9300  Hz  it  is  rated  at 
26  dB  signal  to  noise  ratio.  The  input  signal  range  required  for  a 
normal  level  of  recording  is  specified  at  0.3v  RMS  to  3.0v  RMS  with 
Harmonic  Distortion  rated  at  1%,  tested  at  1  KHz  and  60  ips.  The  input 
impedance  is  greater  than  100  K  ohms.  Each  reproduce  amplifier  has  two 
equalization  plug-in  cards,  with  the  correct  one  being  selected  by  the 
speed  control  switch. 

During  the  playback  operation  the  15  ips  and  1-7/8  ips  speeds  are 
used.  At  15  ips,  the  audio  comments  concerning  amplifier  gains,  maneu¬ 
vering,  calibration  and  general  observations  and  comments  are  monitored 
along  with  a  time  check  of  WWV  compared  to  the  pre-recorded  time  code. 
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The  speed  is  reduced  to  1-7/8  ips  to  give  an  8  to  1  reduction  of 
recorded  frequencies  to  enable  better  reproduction  and  data  acquisition 
by  the  Visicorder.  The  machine  is  operated  manually  and  is  remotely 
controlled  by  the  Tape  Search  Unit. 

The  Systron-Dormer  Model  8130  Time  Code  Reader  is  used  to  translate 
the  recorded  IRIG  B  time  code  (BCD  format)  in  terms  of  hours,  minutes 
and  seconds.  The  incoming  time  code  is  visually  presented  on  a  six 
digit  in-line  nixie  tube  readout  on  the  front,  panel.  The  rack  mounted 
unit  decodes  the  incoming  time  code  in  both  forward  and  reverse  directions 
during  the  data  playback  and  is  used  in  combination  with  the  Model  8140 
Tape  Search  Unit,  thus  providing  a  fully  automatic  tape  search  system, 
with  manual  override  0/ 

The  Systro.i-Donner  Model  8140  Tape  Search  Unit  (TSi)  s  a  combina¬ 
tion  tape  search  and  transport  control  unit  used  with  the  Model  8130 
Time  Code  Reader.  The  TSU  provides  the  remote  control  commands  for  the 
Honeywell  Model  5600  Magnetic  Tape  System  in  order  to  automatically 
search  for,  locate  and  produce  selected  segments  of  test  data. 

The  Model  8140  operates  as  a  digital  limit  comparator  to  compare 
the  decoded  tape  time  from  the  Model  8130  with  selected  start  and  stop 
times  that  are  pre-set  by  the  operator  by  means  of  thumbwheel  switches. 
Upon  detection  of  the  desired  time  comparisons,  the  search  unit  auto¬ 
matically  provides  the  various  control  commands  (contact  closures)  to 
the  tape  transport.  Manual  search  is  also  included  in  the  control. 

^/information  from  Systron-Donner  Data  Products  Division  Manual  for  the 
Model  8130  Time  Code  Reader  for  IRIG  B. 
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In  a  typical  operation  the  operator  sets  in  predetermined  times 
which  are  correlated  with  a  determined  location  relating  to  a  particular 
sample  site  in  the  survey,  i.e.,  the  coring  sites.  The  tape  automatically 
searches  out  the  start  time  and  reads  out  data  to  the  end  of  the  pre¬ 
sented  time  window  at  designated  speeds.  In  order  to  accomplish  this 
feature  a  Variable  Time  Drive  Unit  had  to  be  designed  and  built.  The  MMTC 
modification  to  incorporate  the  Variable  Time  Drive  Unit  (VTDU)  is  used 
to  control  the  operation  of  the  visicorder  paper  drive  motor.  Both  the 
repetition  rate  and  the  duty  cycle  times  can  be  varied  by  the  operator. 
This  unit  was  used  to  reduce  the  amount  of  photographic  paper  used, 
thus  reducing  the  cost  of  playback. 

iypical  sparker  firing  rates  are  .75  second  and  1  second.  The 
visicorder  paper  speed  is  10  inches  per  second.  However,  there  is  an 
8  to  1  reduction  of  tape  speed,  from  15  inches  per  second  to  1-7/8  ips, 
wmch  lengthens  the  firing  rates  to  6  and  8  seconds  respectively.  Paper 
consumption  for  a  10  second  observation  of  10  bursts  amounts  to  600  to 
800  inches  per  site. 

Since  most  of  the  significant  information  could  be  obtained  in  an 
interval  of  .25  second  to  .4  second  the  VTDU  reduced  paper  consumption 
to  a  minimum  of  200  inches  and  to  a  maximum  of  400  inches  per  site. 

The  Honeywell  Accudata  117-06  is  a  rack  mounted  6  channel  direct 
coupled  amplifier  unit  consisting  of  six  separate  solid  state,  wide 
band,  single  ended,  low  gain  amplifiers.  This  instrument  is  used  to 
drive  the  galvanometers  of  the  Model  1508A  Visicorder. 

The  amplifiers  amplify  AC  and  DC  voltages  and  provide  the  propel 
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impedance  matching  and  isolation  between  the  tape  recorder  reproduce 
amplifier  outputs  and  the  visicorder.  Each  channel  has  a  multi-position 
gain  switch  in  steps  of  1,  2,  5  from  .01  to  10  yielding  a  1000:1  range, 
and  a  vemiar  control  that  increases  the  fixed  gain  from  unity  by  a 
factor  of  2.5.  The  maximum  rated  input  is  +_  300  VDC  or  peak  AC.  The 
maximum  rated  output  is  ^  3  VDC  or  peak  AC  at  40  milliamperes . 

The  accuracy  at  DC  is  +_  2%  and  the  gain  or  sensitivity  linearity  is 
+_  0.2%  of  full  scale  output  at  DC.  The  frequency  response  is  flat  within 
+.  5%  from  DC  to  20  KHz.  The  overload  recovery  time  is  less  than  100  micro¬ 
seconds  to  within  +_  1%  of  the  final  value  with  an  overload  10  times  the 
full  scale  input  up  to  a  maximum  of  300  volts.  The  DC  stability  is  rated 
at  less  than  100  microvolt/degree  Centigrade.  The  noise  level  is  loss 
than  +_  2  millivolts  RMS  over  DC  to  20  KHz  for  all  gains. 

The  Honeywell  Model  1508A  Visicorder  (Visual  Recorder)  is  a  direct- 
writing  oscillograph  recorder  with  provisions  for  up  to  12  channels  of 
input.  A  photographic,  permanent  recording  of  the  tape  recorded  analog 
data  is  produced  (figure  3-20f) .  A  mercury  vapor,  ultraviolet  light 
source  produces  the  nominal  0.2  inch  wide  traces  on  the  8  inch  wide 
sensitized  paper  which  is  relatively  insensitive  to  normal  room  lighting. 

The  recorder  includes  a  high-linearity  optical  system  consisting  of 
5  fluid  damped  miniature  galvanometers,  an  accurate  grid  line  system  and 
timing  line  system.  The  recorder  drive  system  can  be  operated  at  speeds 
ranging  from  0.1  to  80  inches  per  second.  Four  reference  timing  line 
intervals  of  .01  to  iO  seconds,  with  a  basic  +.  1.5%  accuracy,  can  be 
selected  on  the  recordings.  The  simultaneous  recording  of  the  grid  lines 


124 


and  galvanometer  traces  at  the  same  point  on  the  paper  eliminates  any  in¬ 
accuracies  due  to  dimensional  changes  in  the  paper  caused  by  humidity 
or  lateral  paper  movement  during  the  recording  process. 

Five,  fluid  damped,  miniature,  type  M-5000  galvanometers  are  used. 
Their  normal  undamped  natural  frequency  is  5000  Hz,  with  a  rated  flat, 

5 %,  frequency  response  of  0  to  3000  Hz,  The  mirrors  are  optically 
coated  to  insure  maximum  efficiency  of  the  reflecting  surface.  The 
maximum  deflection  of  the  AC  operation  is  rated  at  3.5  inches  which 
maintains  the  proper  optical  qualities  and  a  +_  2%  linearity  specification. 
Each  magnet  assembly  is  heated  and  thermostatically  controlled  at  120°  F. 
This  temperature  control  insured  top  performance  and  stability  of  the 
fluid-damped  galvanometers  by  maintaining  a  constant  viscosity  of  the 
damping  fluid. 

A  typical  paper  drive  speed  of  10  inches  per  second  is  used  with  a 
timing  grid  of  0.1  seconds.  Kodak  type  1895  Linagraph  paper  is  used. 

Hie  relative  writing  speed  of  the  paper  is  rated  at  up  to  100,000  inches 
per  second. 

Exact  calibration  of  the  recorded  and  reproduced  information  signals 
is  accomplished  by  recording  calibration  levels  of  1  KHz  sine  wave  signals 
on  each  of  the  active  information  channels  during  the  actual  shipboard 
operation.  The  calibration  levels  are  then  played  back  on  the  shore 
based  equipment  and  any  discrepancy  is  measured,  noted  and  included  in 
all  the  readings  measured  and  computed  from  the  visicorder  photographic 
record.  Therefore,  the  total  error  of  both  the  recording  and  playback 
systems  is  reflected  in  this  one  measured  value  for  each  information 
channel . 
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The  equipment  used  to  record  and  measure  the  calibration  levels  is 
a  Wavetek  Model  111  Voltage  Controlled  Generator  and  a  Tektronix  type 
422  oscilloscope. 

The  various  incoming  information  signals  and  calibration  signals 
are  monitored  with  the  Tektronix  type  422  dual  channel  oscilloscope. 

The  vertical  deflection  accuracy  is  within  the  rated  3%  of  the  indicated 
deflection  factor.  The  range  of  vertical  deflection  is  from  10  milli¬ 
volts  per  division  to  20  volts  per  division,  in  11  calibrated  steps.  The 
frequency  response  of  both  channels  is  flat  from  the  rated  DC  to  15  MHz 
at  the  -3  dB  point.  Single,  dual  trace  alternate,  or  chopped  modes  of 
presentation  are  depicted  for  later  analysis. 

Seismic  Data  Processing. 

All  the  seismic  data  processing  is  accomplished  at  the  Tiburon 
Center.  Visicorder  reproductions  are  analyzed  by  measuring  the  peak 
amplitudes  of  the  transmitted  and  received  pulses.  These  measurements 
are  accomplished  by  hand  and  computed  directly  without  the  aid  of 
digital  computers. 

The  location  of  each  of  the  16  core  sites  are  correlated  for  position 
to  the  seismic  data  by  the  time  reference  marks  on  the  x-y  plots  and  the 
listing  from  the  digital  printer.  The  times  obtained  thusly  are  set  by 
the  operator  on  the  tape  search  unit.  When  the  selected  time  on  the  tape 
is  read  and  monitored,  the  recorder  is  then  set  to  1-7/8  ips  and  the  play¬ 
back  amplifier  adjusted  to  provide  the  proper  recording  levels.  The  tape 
playback  is  automatically  or  manually  cycled  and  by  use  of  the  variable 
time  drive  unit  the  visicorder  drive  is  controlled  to  provide  a  visual 
record  of  the  taped  signals  (figure  3-29f).  Each  tape  is  then  run  to 
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the  recorded  calibration  sections  where  any  signal  level  variations  from 
the  known  calibration  levels  can  be  noted  and  incorporated  in  the  compu¬ 
tation  of  the  measured  signal  levels  (peak  amplitude  of  the  desired  event 
on  the  impulse-response  wave  train).  The  visicorder  gain  levels  are  noted 
on  the  record  (figure  3-29f)  and  the  signal  levels  of  the  desired  events 
are  measured.  The  measurement  is  converted  to  the  corresponding  voltage 
level  for  computation. 

The  steamer  hydrophone  events  are  identified  in  several  ways. 
Generally  the  water  depth  (in  time)  is  sufficient,  or  when  readily  identi¬ 
fied  the  character  of  the  impulse  hydrophone  signal..  Another  method  in¬ 
volves  checking  the  seismic  record  and  timing  the  multiples  of  the  primary 
reflector  (sea  floor  horizon).  One  can  usually  distinguish  the  real 
signal  of  the  surface  reflector  horizon  from  other  reflecting  sources  or 
induced  noise  sources  by  establishing  the  time  differential  between  suc¬ 
cessive  multiple  reflections  and  regressing  to  the  primary  reflector 
event . 

The  playback  wave  forms  fed  from  the  tape  recorder  to  the  Honeywell 
AC  to  DC  amplifier  are  illustrated  in  figures  3-22  and  3-29.  The  typical 
playback  signal  amplitudes  are:  streamer  receiving  hydrophone  signals 
3  to  4vp-p  (figure  3-29a);  the  piggyback  impulse  signals  1.2vp-p  (figure 
3-29b) ;  the  fathometer  receiver  signals  0  to  2vp-p  (figure  3 -29d ) ;  and 
the  fathometer  transmitter  signals  1.5vp-p  (figure  3-29e ) . 

High  Frequency  Inpulse  Signal. 

Trouble  was  encountered  with  both  the  amplifying  and  recording  of 
the  fathometer  received  echo  signals.  The  beam  width  of  the  fathometer 
transducer  is  apparently  too  narrow,  the  fathometer  is  designed  to 
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operate  at  depths  up  to  2460  feet,  to  obtain  an  absolute  amplitude 
measurement  cf  the  signal.  The  primary  problem  encountered  was  main¬ 
taining  a  consistent  return  amplitude  at  any  one  location.  Even  when 
the  boat  was  moored  at  one  fixed  location  to  test  the  consistency,  re¬ 
peatable  amplitudes  were  never  reliable.  For  this  reason,  data  used 
in  the  investigation  is  limited  to  the  broad  band  (0-5  KHz)  sparker 
source  (section  3.2.9). 

3.2.6  Sediment  Sampling  Equipment 

In  order  to  verify  reflectivity  measurements  it  was  necessary  to 
have  a  core  sampler  which  could  rapidly  sample  the  first  meter  of  sur¬ 
face  sediments  and  obtain  relatively  undisturbed  samples.  This  was  re¬ 
quired  for  reliable  engineering  and  mass  physical  property  determinations. 
A  one  meter  core  length  was  decided  upon  because  the  depth  of  penetration 
was  also  within  the  wave  length  of  either  of  the  two  acoustic  sources 
employed.  This  makes  it  possible  to  determine  to  what  degree  changes  in 
the  near  surface  lithology  or  structural  variations  influence  the  re¬ 
flection  coefficient.  Three  samplers  of  the  gravity-type  or  piston- type 
gravity  corer  variety  were  tested  but  found  unsatisfactory.  It  was 
therefore  necessary  to  design,  construct  and  test  a  coring  device  that 
would  provide  the  desired  capability.  It  was  also  decided  that  the 
energy  source  for  penetration  should  not  be  percussion  or  rotary  since 
these  types  tend  to  induce  sediment  distortion.  Therefore, a  steady 
thrust  energy  produced  by  pneumatic  power  was  incorporated  in  the  design. 
Diver  Operated  Corer  -  Design. 

Simplicity  of  operation  and  a  minimal  number  of  mechanical-electri¬ 
cal  features  were  other  criteria  used  for  design  of  the  corer.  A  study 
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of  figure,  3-31  through  3-3-  gives  s.l  f. enplane, ory  details  of  the  lesi.n 

and  hardware  incorporated  in  the  final  product.  Details  covering  the 

laboratory  and  field  test  of  the  prototype  too,  are  in  a  report  by 

Jenkins  and  Takeyama  (1971).  The  field  testing  and  evaluation  of  the 

Pr0t°ty,>e  . . .  conjunction  with  the  sampling  program  for 

the  reflectivity  experiment. 

Operation. 

AH  work  associated  with  the  coring  program  is  done  fro*  the 

DOODlEBUG,  a  40-foot  converted  launch  (sect  3  2  71  a 

^sect.  i. 2.7).  A  crane  and  winch, 

having  a  load  capacity  of  1,000  pounds,  were  used  to  handle  the  corer. 
scuba  divers  are  usdd  during  the  actual  sampling.  three  divers  are 
normally  required;  two  operate  the  device  under  water  and  the  third 
remains  on  ,h.  surface  as  standby  and  relief.  Vaives,  located  on 

. . .  *  th'  Tb.se  can  be  regulated  by  the  diver 

to  enable  a  flow  of  air  pressure  to  drive  a  piston  that  forces  the 

Shelby  sampling  tube  slowly  l„,0  . . floor.  Th.  force  of  the  steady 

thrust  Is  controlled  by  regulating  the  pressure,  when  full  penetration 
is  observed  a  three-way  valve  is  reversed  and  three  independent  hydraulic 
pistons  lift  ,h.  main  carriage  away  from  th.  base  pl«.,  thus  retrieving 
the  core  barrel  from  the  sediment.  A  diver  Immediately  caps  th.  ,„w.r 
.nd  of  th.  barrel  to  preclude  core  loss.  The  cor.  uni,  Is  complete,, 
sealed  while  under  water  so  that  water  content  remains  as  „  „,s  situ 

The  removal  of  ,h.  cor.  barrel  and  replace™,  of  a  second  cor.  barrel  Is 
done  by  the  divers  under  water.  . . .  apd  ^  preparat,on 

for  the  ,i«t  sampling  operation  more  efficient,  particularly  when  a 
vessel  has  a  very  small  deck  area. 
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Figure  3-32.  Diver  operated  oorer  (extended  position) . 
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Figure  5-33.  Core  retainer  modified  Shelby  tube. 
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A  complete  cycle  of  sampling, excluding  ship  positioning  and  2-point 
mooring  time, takes  approximately  10-20  minutes  in  water  40-60  feet  deep. 
The  corer  is  capable  of  sampling  both  semi-consolidated  (clays)  as  well 
as  unconsolidated  sediments. 

3.2.7  Research  Vessel 

All  the  fiflld  operations  mentioned  in  this  experiment  were  conducted 
from  the  R/V  Doodlebug,  a  40-foot  converted  launch  (figure  3-35).  The 
vessel  is  equipped  with  standard  shipboard  equipment  consisting  of  a 
compass,  single  side  band  radio,  radar,  fathometer,  a  4  KW  diesel-electric 
generator  to  supply  110/220  VAC  power  to  '’he  recorders  as  well  as  a 
separate  7.5  «CW,  110/220  VAC  diesel-electric  generator  to  operate  the 
sparker  seismic  system.  The  vessel  is  powered  by  two  197  hp  diesel 
engines  with  twin  screws  for  easy  maneuverability.  Cruising  speed  is 
about  9  knots,  although  normally  survey  speeds  average  4  to  5  knots.  A 
1/2-ton  portable  hoist  with  600  feet  of  3/16"  stainless  steel  cable  is 
mounted  in  the  aft  section  for  handling  loads  up  to  1000  pounds  or  towing 
vehicles  containing  instruments. 

Precision  Navigation  and  Positioning  System. 

The  precision  radio  location  system  is  manufactured  by  Hastings- 
Raydist,  Inc.-  (subsidiary  of  Teledyne,  Inc.).  It  is  used  in  referencing 
ship's  position  to  geodetic  control  and  for  automatic,  real-time  tracking 
during  survey  operations.  Incorporated  in  the  system  is  an  electronic 
subsystem  which  enables  real-time,  on-line  digital  and  graphical  con¬ 
version  of  lane  identification  coordinates  derived  from  the  Raydist 
range-range  mode  system  to  x-y  grid  coordinates.  The  subsystem  was  de¬ 
veloped  by  the  Marine  Minerals  Technology  Center.  Diagrams  and  hardware 
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components  of  the  total  system  are  shown  in  figures  3-36  and  3-37, 

The  ship's  position  and  track  are  accurately  located  to  within  +  10  feet 
based  on  daily  repeatability  checks  of  the  calibrate  value  (geodetic 
control  point).  Details  concerning  the  system  and  electronic  subsystem 
may  be  found  in  a  report  by  Barnes  and  Newman  (1971). 

3.2.8  Field  Operations 

The  field  test  site  selected  for  this  investigation  was  originally 
intended  for  Monterey  Bay,  California;  however  time  and  circumstances 
prevented  accomplishing  this  objective.  The  investigation  reported 
herein  is  confined  to  a  site  in  San  Francisco  Bay  which  had  been  intended 
for  use  as  a  preliminary  site  for  shakedown  tests  of  the  systems  and  pro¬ 
cedures  prior  to  the  Monterey  Bay  study.  Although  a  thorough  and  de¬ 
tailed  study  as  programmed  for  the  Monterey  Bay  work  was  not  conducted  in 
the  San  Francisco  Bay  experiment,  it  is  felt  the  data  obtained  is  signifi¬ 
cant  and  representative  of  the  goals  desired  and  worthwhile  documenting  in 
this  report  (section  3.2.9). 

Target  Site. 

An  area  immediately  south  and  west  of  Red  Rock  in  the  northern 
reaches  of  San  Francisco  Bay  and  south  of  the  San  Rafael -Richmond  P' Adge 
was  selected  for  the  test  site  (figure  3-38).  The  bottom  sediments  in 
the  survey  area  are  of  several  types.  Mud,  silt,  sand,  rock  and  clays 
indurated  with  shells  comprise  some  of  the  types  of  materials  dredged 
from  the  bottom  in  earlier  work  by  the  National  Ocean  Survey  (formerly 
U.  S.  Coast  and  Geodetic  Survey  --  refer  to  chart  USCGGS  No.  5532). 

The  wide  variability  of  sediment/rock  type,  shallow  water,  estuarine 
environment,  and  logistics  were  factors  influencing  the  selection  of 
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Figure  3-36 .  Block  diagram  basic  navigation 
system  and  subsystem. 


'  'Ijish'n  tir* 


Figure  3-38.  San  Francisco  Bay  test  site  (proposed  program). 


the  area  and  determining  the  core  sites  (see  Appendix  -  Sect.  4.5) 

A  program  was  drawito  apply  a  dense  grid  pattern  of  600  foot  line 
spacing  over  an  area  of  1-1/2  square  miles  (figure  3-38).  Core  locations 
were  designated  at  16  sites,  each  located  approximately  900  feet  apart. 
The  actual  program  proposed  was  completed  almost  in  total  and  as  pro¬ 
grammed.  Two  separate  seismic  surveys  were  made  at  different  times  in 
August  1971  and  September  1971.  Figure  3-39  reflects  the  latter  survey 
covering  the  actual  work.  Most  of  the  reflectivity  data  was  derived 
from  the  second  survey.  Seismic  production  amounted  to  about  26.5  miles 
of  continuous  seismic  profiling.  Sixteen  core  locations  were  occupied 
and  at  least  one  core  was  taken  at  each  site.  At  several  core  sites 
two  or  more  cores  were  taken.  Core  retention  averaged  about  90%  for  all 
sites.  Water  depth  ranged  from  26  to  41  feet  (uncorrected),  averaging 
about  33  feet.  The  bottom  topography  is  generally  uniform  with  only 
minor  changes  occurring  in  the  western  edge  of  the  area  in  the  vicinity 
of  the  Richmond  Channel  (figure  3-38). 

Navigation  Positioning. 

The  Raydist  Model  DR-S  High  Accuracy  Short -to-Medium  Range  Radio¬ 
location  system  was  used  to  provide  precise  positioning  and  guidance  for 
the  vessel  (section  3, 2. 7.1).  Two  shore-based  stations  are  required  for 
the  range-range  system.  The  north  station,  designated  "Red",  is  located 
at  the  University  of  California  Field  Station  in  Richmond  and  the  south 
station,  designated  "Green",  is  located  at  the  Treasure  Island  Naval  Base 
in  the  center  of  San  Francisco  Bay  mid-span  of  the  Oakland-San  Francisco 
Bay  Bridge. 
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Figure  3-39 


San  Francisco  Bay  test  site 


actual  survey. 


Pre-p.ots  of  profile  track  lines  are  drawn  on  transparent  rouseable 
overlays  made  from  conventional  maps  using  the  California  Grid  system 
(Lambert  Coordinates)  as  a  reference.  He  boat  operate.-  uses  the  x-y 
Plotter  with  preplotted  track  to  guide  hi.  onto  the  desired  course  and 
to  pre-designated  core  sites.  The  scale  of  the  overlays  is  large 
enough  to  provide  high  resolution  plotting  capability  (,  i„ch  .  1000  feet) 
A  Hewlett-Packard  Model  9100A  programmable  calculator  automatically  con¬ 
verts  the  range-range  "red"  and  "green"  lane  counts  that  reflect  distance 
from  each  of  the  respective  shore  stations  to  the  vessel  (each  lane  equals 
148.2282  feet)  to  x  and  y  Lambert  grid  coordinates  (California  State  Grid 
Coordinate  System).  The  Hewlett-Packard  Model  912SA  digital  pi„tter  is 
fed  these  coordinates  approximately  every  1.3  seconds  to  provide  a  real¬ 
time  automatic  plotting  navigation  positioning  system.  A  digital  printer 
.s  commanded  through  the  program  to  print  out  on-line  the  "red"  and  "green" 
U"e  counts,  x  and  y  coordinates  and  the  time  of  print  out.  Here  is  also 
a  provision  for  digital  depth  data  to  be  printed  if  it  is  required.  Green¬ 
wich  Mean  Time  is  used  as  the  reference  for  all  events  recorded  in  the 
navigation-positioning  operation,  and  for  ancillary  recorders  associated 
with  the  acoustic  and  magnetic  tape  recording  systems. 

He  average  position  error  for  the  runs  made  on  a  daily  basis  was 
•  06  lane  on  the  "red"  count  (8.9  feet)  and  .08  lanes  on  the  "green"  count 
(11.9  feet).  He  Raydis.  is  a  non  line-of-sight  system  using  ground  wave 
transmission.  Errors  are  felt  to  be  caused  by  changes  in  velocity  of  wave 
propagation,  not  positioning  the  boat  in  the  exact  position  when  it  is 

originally  calibrated  and  recalibrated,  or  insufficient  warmup  time  prior 
to  commencing  operations. 
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The  subsyttem  interface  also  contains  a  digital  clock  with  5  select¬ 
able  time  outputs  ranging  from  1  minute  to  30  minutes.  This  timer  controls 
a  marker  interface  which  simultaneously  marks  the  records  of  the  fathom¬ 
eter,  seismic  recorder,  Raydist  strip  chart  recorder,  and  controls  the 
Raydist  digital  printer  (auxiliary unit  to  the  HP  printer).  The  vessel's 
location  in  lane  counts  is  sequentially  assigned  station  numbers  that  are 
used  for  annotation  in  addition  to  the  time  reference.  The  operator 
periodically  annotates  the  records  with  the  station  number  and  time.  The 
marker  interface  can  also  be  manually  operated  to  index  significant  events, 
such  as  occur  during  maneuvering,  amplifier  gain  changes,  calibration 
points,  coring  events,  and  so  forth. 

Survey  and  Sanpling  Operations. 

The  total  complement  of  crew  and  project  staff  consists  of  a  boat 
operator,  an  electronics  technician  and  a  geophysicist  during  seismic 
operations.  The  sampling  work  requires  three  divers  and  a  geologist  in 
addition  to  the  boat  operator  and  electronics  technician. 

Once  the  navigation  system  is  calibrated,  at  a  preselected  geodetic 
control  point,  the  vessel  proceeds  to  the  test  area.  The  acoustic  source 
and  receiver  units  are  deployed  and  towed  in  a  side-by-side  configuration 
approximately  65  feet  astern  of  the  vessel.  Vessel  speed  during  survey 
operations  is  maintained  at  about  4  -  S  knots.  The  seism:! c  source  is 
electronically  detonated  once  every  second  and  the  received  signal  is 
recorded  both  in  graphic  form  (variable  density)  and  onto  magnetic  tape 
(section  3. 2.5. 3).  Voice  annotation  on  the  tape,  and  written  logs  are 
maintained  throughout  the  operations. 
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During  sampling  operations,  the  vessel  is  anchored  and  position 
movement  is  monitored  on  the  x-y  plotter  readout  (scale  1"  =  50'). 
coring  is  accomplished  at  slack  tide  to  allow  safe  operations  for  the 
divers.  Core  production  was  limited  to  two  per  day  mainly  because  of 
the  strong  currents  produced  by  the  changing  tide.  Details  of  the  coring 
operation  are  noted  in  section  3. 2. 6. 2. 

3.2.9  Data  Analysis 
Core  Sample  Analysis. 

The  object  of  the  core  sample  aniysis  phase  of  this  investigation  is 
to  determine  the  character  of  the  sediments  in  the  area  where  the  acoustic 
measurements  are  made.  The  analyses  are  directed  toward  obtaining  grain- 
size  and  mass  characteristics  of  the  sediments.  These  measurements  are 
used  directly  for  correlation  with  the  acoustic  reflectivity  measurements. 
Due  to  the  shortage  of  time  and  personnel,  only  cursory  mineralogical 
examinations  have  been  made  of  the  samples  at  the  time  of  this  report. 

The  samples  were  obtained  with  the  WfTC-developed,  diver-operated, 
bottom-sitting  coring  drill,  which  is  capable  of  obtaining  a  core  70 
centimeters  long  and  7.9  centimeters  in  diameter  (sections  3.2.6  and 
3. 2. 8. 3).  Most  of  the  cores  obtained  during  this  investigation  were 
between  60  and  65  centimeters  in  length.  When  the  sample  is  taken,  the 
core-barrels  are  sealed  under  water,  brought  to  the  surface,  and  stored 
in  a  vertical  position  during  transit  to  the  Tiburon  laboratory.  Usually, 
the  laboratory  processing  begins  within  2  to  4  hours  after  the  sample  is 
ta!.en  from  the  floor  of  San  Francisco  Bay. 

Upon  arrival  in  the  laboratory,  the  exterior  of  the  core-barrel  is 
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cleaned  and  wiped  dry,  and  the  whole  assembly  (core-barrels,  seals  and 
contained  sample)  weighed.  The  seal  is  then  removed  and  the  excess 
water  carefully  siphoned  off  and  it3  weight  recorded.  The  depth  to  the 
solid  material  is  then  measured  and  recorded,  and  the  volume,  wet  weight, 
and  wet  density  of  the  sample  calculated.  The  core-barrel  is  then  placed 
in  the  horizontal  core-extruder  and  extrusion  of  the  core  starts. 

Samples  for  water  content  determinations  are  obtained  by  taking  six 
slices  1  cm  thick  from  the  core  as  it  is  extruded  from  the  barrel.  These 
slices  are  caught  in  pre-weighed  glass  Petri  dishes,  and  each  set  of  six 
samples  is  weighed  immediately,  placed  in  a  calibrated  drying  oven,  and 
dried  to  constant  weight.  The  process  of  obtaining  the  slices,  weighing, 
and  placing  the  samples  in  the  drying  oven  consumes  approximately  3 
minutes  per  set  of  6  samples.  After  the  material  for  the  water  content 
determinations  is  obtained,  Torvan^and  Penetrometer*^ measurements  are 
made  on  the  face  of  the  core  remaining  in  the  core-barrel.  The  core  is 
extruded  another  14  centimeters,  cut  off,  and  stored  in  plastic  containers 
for  future  mineralogical  analysis.  This  routine  is  repeated  until  no 
further  core  remains  in  the  core-barrel. 

Specific  gravity  (dry  density)  determinations  are  made  on  the 
water-content  samples  after  they  have  been  dried  and  weighed.  A  50  -  60 
gram  split  of  the  dried  material  is  taken,  crushed,  and  ground  to  approxi¬ 
mately  minus  100-mesh  size  and  the  specific  gravity  is  obtained  with  a 


‘The  Torvane  and  Penetrometer  data  are  not  included  in  this  report 


because  of  suspected  instrument  malfunction. 


Beckman  Model  i>30  Air  Comparison  Pycnometer  for  each  of  three  levels  in 
the  core;  0-6  cm,  20-26  cm,  and  40-46  cm  with  0  corresponding  to  the 
surface  of  the  bay  floor.  The  average  of  the  three  determinations  is 
used  as  the  dry  density  for  the  core  sample. 

‘Hie  remainder  of  the  water-content  material  is  placed  in  a  200-mesh 
sieve  and  the  minus  200-mesh  material  is  washed  from  the  sample  into  a 
3-gallon  plastic  bucket,  allowed  to  settle,  and  the  clear  water  carefully 
decanted.  The  slurry  remaining  in  the  bucket  is  washed  into  a  drying  pan 
with  a  minimum  amount  of  water  and  dried.  The  dried  minus  200-mesh 
material  is  weighed,  and  a  35  gram  split  taken  for  sub-sieve  analysis. 
Sub-sieve  analysis  is  accomplished  by  the  hydrometer  method  using  1000  ml 
graduate  cylinders  in  a  constant  temperature  water  bath.  Meanwhile,  che 
plus  200-mesh  material  is  dried  ar.-*  sieved  through  selected  sieves  using 
a  mechanical  sieve  shaker.  A  10-  minute  shaking  cycle  is  used  for  the 
sieve  analysis.  A  per  cent  finer  cumulative  curve  is  constructed  from 
the  sieving  and  hydrometer  analyses  and  the  median  grain  size,  25  per 
cent  quartile,  and  the  75  per  cent  quartile  picked  off  the  curve.  The 
arithmetic  quartile  deviation,  the  geometric  quartile  deviation  (sorting 
coefficient  of  Trask)  and  the  Log10  Sorting  Coefficient  are  then  calcu¬ 
lated.  A  tabulation  of  the  grain  size  and  mass  characteristics  of  all 
the  sediment  samples  is  presented  in  the  Appendix  (section  4.6). 
Correlation  of  Mass  Physical  Properties  to  Reflectivity. 

The  grain  size  -nd  mass  properties  data  are  correlated  with  the 
acoustic  data  using  a  linear  regression  program  on  a  Hewlett-Packard 
Model  U100A  Programmable  Desk  Top  Calculator,  ihc  program  calculates 
the  equation  of  the  straight  line  of  best  fit  of  a  set  of  data  points. 


145 


The  best  fit  is  determined  by  minimizing  the  sum  of  the  squares  of  the 
deviations  of  the  data  points  from  the  line.  The  program  also  calculates 
a  correlation  coefficient,  r,  an  indication  of  goodness  of  fit.  If  r  =  0, 
there  is  no  correlation.  If  r  *  ♦_  1,  there  is  perfect  correlation  or  a 
perfect  fit.  Correlations  were  run,  where  data  available,  on  the  0-6  cm 
fraction,  20-26  cm  fraction,  40-46  cm  fraction,  and  also  on  the  average 
value  for  the  whole  core  of  each  of  the  grain  size  and  mass  characteristics 
obtained.  The  results  of  these  correlations  are  shown  in  Tables  3-3  and 
3-4,  as  well  as  the  standard  error  per  correlation. 

3.2.10  Results  and  Conclusions 

The  magnitude  of  the  correlation  coefficients  presented  in  Tables 
3-3  and  3-4  corroborate  Breslau's  (1967)  conclusion  that  a  significant 
correlation  exists  between  bottom  loss  and  certain  physical  characteristics 
of  the  sediments.  Most  notable  are  the  relationships  between  bottom  loss 
and  the  sediment's  porosity,  grain  size,  and  fine  material  (silt-clay) 
content.  A  comparison  of  the  correlative  values  obtained  in  this  investi¬ 
gation  and  those  obtained  by  Breslau's  (1967)  much  more  extensive  investi¬ 
gation  is  presented  in  Table  3-5.  It  should  be  noted  that  Breslau's  data 
were  obtained,  for  the  most  part,  in  the  offshore  marine  environment  where 
relatively  uniform  sediments  prevail  over  large  areas,  whereas  this  in¬ 
vestigation  was  conducted  in  essentially  an  estuarine  environment  wherein 
the  character  of  the  bottom  sediments  is  extremely  variable  over  very 
short  distances.  The  comparability  of  the  two  sets  of  data  is  therefore 
rather  remarkable  and  seemingly  indicates  a  high  degree  of  validity  for 
the  technique.  Also,  the  improvement  must,  in  part,  be  attributed  to 
the  use  of  precision  radiolocation  positioning  and  navigation  equipment 


146 


10 

Table  3-3.  Correlation  Between  Bottom  Loss - 
and  Sediment  Grain-Size 


/ 


Sediment  Grain  Size 
(core  interval) 


Correlation 
Coefficient  (r) 


Standard  Error— ^ 
dB/Unit 


Median  Grain  Size  (0-6cm) 

.87 

♦  .325 

Median  Grain  Size  (20- 26cm) 

.83 

♦  .367 

Mediam  Grain  Size  (40-46cm) 

.91 

♦  .276 

Median  Grain  Size  (av  for  ea  core) 

.92 

1  -227 

%  (Silt  6  Clay)  (0-6cm) 

.93 

♦  .252 

%  (Silt  &  Clay)  (20-2Ccm) 

.75 

♦  .450 

%  (Silt  6  Clay)  (40-46cm) 

.89 

+  .308 

%  (Silt  8  Clay)  (av  for  ea  core) 

.89 

♦  .305 

Shell  content  (0-6cm) 

.004 

♦  .643 

Shell  content  (20-26cm) 

.27 

♦  .692 

Shell  content  (40-46cm) 

.13 

♦  .609 

Shell  content  (av  for  ea  core) 

.69 

jf  .4651 

Sorting  Coefficient  (0-6cm) 

.89 

♦  .291 

Sorting  Coefficient  (20-26cm) 

.78 

+_  .395 

Sorting  Coefficient  (40-46cm) 

.83 

♦  .356 

Sorting  Coefficient  (av  for  ea  core) 

.92 

♦  .252 

-  Bottom  loss  calculations  are  based  on  peak  pressure  amplitude  measure¬ 
ments  and  Rayleigh  Reflection  Coefficient  determinations. 

11/  The  standard  error  of  estimate  (S)  is  the  standard  deviation  about  a 
line  of  average  relationship.  Assuming  a  normal  distribution,  68  per  cent 
of  all  cases  will  lie  within  a  range  of  *_  S,  95  per  cent  will  fal’  within 
♦  2S,  and  99.7  per  cent  will  fall  within  ♦  3S.  In  this  instance,  the  re¬ 
lationship  of  bottom  loss  to  the  other  variable  is  of  interest.  The  values 
of  S  are  therefore  in  terms  of  decibels  (db)  per  unit  of  paired  variable. 
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and  developing  a  sampler  that  can  obtain  relatively  undisturbed  core 
samples . 

Figures  3-40  and  3-41  provide  graphic  illustration  of  the  practical 
value  of  the  investigation.  The  3-component  diagram  (figure  3-40)  shows 
a  marked  propensity  toward  finer  grained  material  in  the  San  Francisco 
Bay  test  site.  It  is  significant  to  also  note  the  coincidence  of  the 
data  along  the  western  edge  to  the  lafge  area  designated  for  dumping  of 
spoils  indicated  in  the  nautical  chart  (USCfjGS  Chart  5532).  The  high 
degree  of  fines  may  be  associated  with  the  dump  area  (figure  3-41). 

In  producing  maps  from  acoustic  measurements  such  as  represented  by 
figure  3-41,  it  is  important  to  realize  how  meaningful  the  distribution 
factors  are.  The  lithologic  designations  (i.e.,  sand,  clay,  silt,  etc.) 
illustrated  on  the  map  would  appear  at  first  sight  to  reflect  only  gross 
classification,  without  regard  to  their  respective  mass  physical  properties. 
However,  the  addition  of  bottom  loss  (dB)  characterizes  much  more  than 
lithology.  Within  the  zone  of  influence  there  is  an  added  dimension  that 
can  provide  the  practical-oriented  engineer  meaningful  data  regarding 
soil  and  rock  properties.  The  combination  can  yield  valuable  information 
leading  up  to  the  prediction  of  parameters  such  as  stability,  bearing 
strengths,  compressibility,  etc. 

Last,  but  not  least,  this  knowledge  can  provide  the  ability  to  under¬ 
stand  how  disturbance  of  these  materials  (i.e.,  through  turbidity,  trans¬ 
portation  of  fines  in  solution)  will  affect  the  environment  and  influence 
the  ecologic  balance. 
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Table  3-4.  Correlation  Between  Bottom  Loss 
and  Sediment  Mass  Physical  Properties 


12/ 


Correlation 
Coefficient  (r) 


Standard  — 
Error  of 
Estimate 


Dry  Density 

.03 

+ 

.679 

Wet  Density 

Water  Content 
Water  Content 
Water  Content 
Water  Content 

(0-6cm) 

(20- 26cm) 

(40 -46 cm) 

(average  for  each 

.70 
.90 
.84 
.86 
core).  88 

+ 

+ 

+ 

+ 

+ 

4.166 

.295 

.374 

.352 

.328 

Void  Ratio 

.68 

+ 

.452 

Porosity 

.73 

+ 

.455 

Table  3-5.  Comparison  of  Correlative  Values  of  MMTC  vs  Breslau  (1967) 

13 

Standard  Error 


Correlation 

MMTC 

Coefficient 

Breslau 

dB/Unit 

MMTC 

Breslau 

Porosity 

.73 

.70 

+  .455 

+  .039 

Grain  Size 

.92 

.646 

+  .227 

+  .292 

Sorting 

.92 

.343 

+_  .252 

+  .660 

Fine  Material 
(Silt  8  Clay) 

.89 

.786 

+  .305 

1  -015 

!*/  Bottom  loss  calculations  are  based  on  peak  pressure  amplitude  measure 
ments  and  Rayleigh  Reflection  Coefficient  determinations. 


'3/  The  standard  error  of  estimate  (S)  is  the  standard  deviatio^about^ 
line  of  average  relationship.  Assuming  a  normal  distribution,  68  per  cen 
of  all  cases  will  lie  within  a  range  of  +  S,  95  per  cent  will  fall 
+  2S  and  99.7  p<rcent  will  fall  within  +  3S.  In  this  instance,  the  re- 
Tationship  of  bottom  loss  to  the  other  variable  is  of  interest  The  values 
of  S  are  therefore  in  terms  of  decibels  (db)  per  unit  of  paired  variable. 
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APPENDIX 


A- I  Data  for  Resistivity  vs  Porosity  Measurements 


Water  _ 

Sample  Conductivity  Porosity  Resistivity  Conductivity  Formation 
I.D.  mho/m  %  ohm-m  mho/m  Factor 


35x48  4.70  42.8 

sand  41.2 

39.8 

38.1 


48x65  4.40  42.7 

sand  39.7 

37.3 


0.757  1.32  3.56 

.782  1.28  3.68 

.807  1.24  3.79 

.836  1.20  3.94 


.806  1.24  3.55 
.900  1.11  3.96 
.945  1.06  4.33 


100x150  4.70 

sand 


46.2 

.674 

43.7 

.713 

42.7 

.723 

40.8 

.802 

39.3 

.812 

1.48 

3.12 

1.41 

3.35 

1.38 

3.40 

1.25 

3.77 

1.23 

3.82 

4x6  4.40 

dolomite 


49.6 

.683 

47.2 

.718 

45.0 

.757 

40.0 

.811 

1.46 

3.00 

1.39 

3.16 

1.32 

3.33 

1.23 

3.57 

9x10  4.40 

dolomite 


49.3 

47.0 

44.0 


.732 

1.37 

3.22 

.786 

1.27 

3.46 

.846 

1.18 

3.86 

14x20  4.25 

dolomite 


51.0 

.743 

49.3 

.786 

47.4 

.832 

1.35 

1.27 

1.20 


3.17 

3.34 

3.53 


65x103  4.65  55.3 

dolomite  54.2 

52.8 

50.8 
49.3 


.584 

1.71 

2.72 

.594 

1.68 

2.83 

.634 

1.58 

2.97 

.673 

1.49 

3.12 

.678 

1.47 

3.15 
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Sample 

I.U. 

Water 

Conductivity 

mho/m 

Porosity 

% 

Resistivity 

ohm-m 

Conductivity 

mho/m 

Formation 

Factor 

35x48 

sand 

3.40 

42.5 

40.7 

39.0 

36.5 

36.1 

1,00 

1.04 

1.12 

1.19 

1.20 

1.00 

0.980 

0.893 

0.840 

0.833 

3.40 

3.55 

3.79 

4.05 

4.08 

48x65 

sand 

3.43 

43.9 

42.2 

40.3 

38.3 

0.965 

1.01 

1.07 

1.14 

1.04 

0.990 

0.934 

0.877 

3.30 

3 .  i!  7 

3.68 

3.89 

100x150 

sand 

3.47 

44.7 

43.1 

41.3 

38.8 

0.955 

0.990 

1.04 

1.11 

1.05 

1.01 

0.962 

0.901 

3.32 

3.44 

3.63 

3.85 

4x6  3.48 

dolomite 


49.0 

0.842 

47.4 

,881 

46.2 

.895 

44.4 

.955 

1.19 

2.93 

1.135 

3.07 

1.12 

3.13 

1  .  U5 

3.23 

9x10  3.48 

dolomite 


50.3 

.88] 

48.5 

.  9  j5 

47.3 

.970 

45.6 

1.02 

44.3 

1.05 

1.135 

3.07 

1.07 

3.26 

1.03 

3.38 

0.980 

3.57 

.953 

3,68 

14x20  3.48 

dolomite 


51.6  .891 

50.0  .945 

48.4  1.01 

46.1  1.06 


1 • 12  3.11 

1-06  5.30 

0-990  3.51 

•943  3.71 


65x100  3.48 

dolomite 


54.5 

.827 

52.0 

.852 

50.8 

.881 

48.3 

.955 

1 • 21  2.88 

1  •  17  2.97 

1-135  3.07 

1  •  05  3.33 
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Water 


Sample 

I.D. 

Conductivity 

mho/m 

Porosity 

% 

Resistivity 

ohm-m 

Conductivity 

mho/m 

Formation 

Factor 

35x48 

1.67 

42.7 

2.05 

0.488 

3.43 

sand 

40.6 

2.16 

.463 

3.56 

38.3 

2  32 

.432 

3.87 

36.3 

2.46 

.407 

4.12 

48x65 

1.67 

44.6 

1.93 

.518 

3.23 

sand 

42.6 

2.08 

.482 

3.46 

40.7 

2.18 

.458 

3.64 

37.4 

2.42 

.413 

4.05 

100x150 

sand 

1.67 

46.3 

43.4 

42.2 

39.2 

1.79 

1.91 

1.98 

2.28 

.559 

.524 

.505 

.438 

2.99 

3.20 

3.30 

3.81 

4X6 

1.67 

50.0 

1.66 

.603 

2.77 

dolomite 

48.5 

1.72 

.582 

2.87 

47.0 

1.78 

.562 

2.97 

43.5 

1.98 

.505 

3.30 

9x10 

dolomite 

1.67 

50.3 

48.7 

47.2 

43.7 

1.73 

1.84 

1.95 

2.15 

.578 

.543 

.513 

.465 

2.89 

3.07 

3.24 

3.59 

4x20 

1.74 

53.3 

1.63 

.613 

2.84 

dol  >mite 

51.9 

1.69 

.592 

2.93 

50.3 

1.78 

.562 

3.10 

48.7 

1.88 

.533 

3.28 

46.5 

1.98 

.505 

3.45 

65x100 

1.76 

55.9 

1.43 

.700 

2.61 

dolomite 

54.1 

1.53 

.654 

2.69 

53.2 

1.55 

.645 

2.73 

51.8 

1.62 

.617 

2.84 

50-0 

1.73 

.578 

3.03 
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A-2 


Calculation  of  the  Potential  Field  Due  to  a  Poin 
Source  of  Current  over  a  Two- layered  Sea  FI 


A-2- I  Introduction 

In  order  to  develop  curves  such  as  those  shown  in  figures  2-14 
and  2-15,  the  potentials  about  a  point  source  of  current  on  the  sea 
floor  must  be  known.  Terekhin  (1962)  hes  derived  e  general  method  for 
calculating  the  potentials  over  a  bottom  with  any  number  of  layers.  an< 
yan  (1956)  has  developed  curves  for  a  limited  nunter  of  cases, 
harke  (1965)  has  done  a  more  detailed  study  of  a  Wenner  array  on  the 
sea  floor  over  a  single  submarine  layer.  |„  this  section,  equations 
for  a  point  source  of  current  over  a  one- or  two- layered  bottom  wi  1 1  be 

developed  and  a  computer  program  for  calculating  these  potentials  will 
be  presented. 

The  necessary  relationships  between  current  density,  conductivity, 
charge  density,  potential,  and  time  are: 


Ohm's  law:  3  =  at  *  £/ p 

Continuity  equation:  s  -3p/3t 

Electric  field  Intensity:  £  *  -$V 

where  iJ  Is  current  density 

a  Is  electrical  conductivity 
p  Is  electrical  resistivity 

^  is  the  gradient  operator, 

(In  rectangular  coordinates) 
t  Is  time 

V  is  potential 


(4.1) 

(4.2) 

(4.3) 
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Combining  these,  and  assuming  the  medium  to  be  homogeneous  and  isotropic, 
g  ives 


a72V  =  3p/3t 

or 

=  Ip 

where 

I  is  the  current 

(4.4) 


2  .  /  3  2 

V  is  the  Lapiacian  operator,  I 

(in  rectangular  coordinates) 

if  no  current  is  being  generated  at  the  point, 
V2V  =  0 


3y‘ 


3Z  1 


■) 


(4.5) 


The  potential  V  at  any  point  in  a  region  in  which  current  Is 
flowing  is  given  by  (4.4),  and,  if  no  current  Is  being  generated  at 
the  point,  the  potential  is  given  by  (4.5).  For  the  special  case  of 
a  whole  space  containing  a  point  source  of  current,  the  solution  to 
(4.4)  is 

V  -  I  p/4ttR  (4.6) 

where  R  is  the  distance  from  the  current  source. 


•  63 


A  -  2  -  2  General  Solution 


In  cylindrical  coordinates  (4.5)  may  be  expressed  av 

V>V(r ,e.z)  .0*IfitfiV*Ir|iV.o 


In  our  proble,  we  will  have  no  0  dependence  so  (4.7)  may  be 
simpl If  led  to 


32V  +  1  J.  +  32V  , 


dr 


r  r  dz  ‘ 


(<♦.8) 


By  separating  variables,  the  solution  to  (4.8)  Is  shown  to  be 


V(r,z)  «J  [A(A)eXz  +  B(x)e"XzJ  Jo(Ar)dA 

w^,er®  A  and  B  are  constants 

X  Is  a  dummy  variable 

Jo  Is  the  Bessel  function  of  order  zero 
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A- 2 - 3  Particular  Solution 


(4.9) 


The  two-layer  case  analyzed  is  shown  In  figure  4-1.  In  any 
layer  not  containing  the  current  source,  the  potential  Is  given  by: 

vi(r.z)  *r[A1(X)eXz  ♦  Bi(X)e‘Xz]  Jo(Ar)dA 

where  l  is  2,  I ,  or  -I.  In  layer  zero,  the  particular  solution 
(4.6)  must  be  added: 

V.  ■  ♦  Wao(X)«X2  +  B.Ule'*2]  Jo(Xr)dX  , 

where  R  » j~ 


Letting 


rZ  +  ZZ 
!P, 

%  * 


the  complete  solution  to  the  problem  may  be  expressed  by  the 
following  set  of  equations: 

V2  «(  [a2 (A)eXz  +  B2 (A )e“Xz]  Jo ( A r )d A 

■'O 

Vi  [A,(A)eXz  +  B  i  ( A )  e_Xz  ]  J0(Ar)dA 

Vo  *  ip*  |ao  (A  )eXz  +  Bo(A)e  Xzj  Jo(Ar)dA 

V-i  ■  j  [a. i ( A)eXz  +  3.I(A)e’Xzj  J0(Ar)dA 

The  boundary  conditions  at  each  interface  are: 


V1  *  V1-, 


1  3V1 


1  3V1-. 


1  *  2,  1,  0 


and 


dz 


V2— 0 

Vi— 0 


Pf-i 


as 

as 


z— ►-  00 


(4.10) 


(4.11) 

(4.12) 

(4.13) 

(4.14) 
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(4.15) 

(4.16) 

(4.17) 

(4.18) 


Crin'  ^  West  <l965>  «Wtlon,  (4.  ,5)  ,„d  (ll.16)  „„ 

pages  404-1*05,  and  Van  Nostrand  and  Cook  (1466)  treat  the  n  hi 

v i yoo)  treat  the  problem  outlined 

above  for  the  land  case  In  detail. 


By  using  the  Weber-L Ipsch I tz  identity, 

i  C 00 

IT  *  ]o  e'XZ  Jo(Ar)dA  ;  z  >  0 

1  f  °° 

R  s  Jo  eU  J*(Ar)dX  ;  z  <  0 

def In ing 


(4.19) 

(*♦.20) 


koi  =  (Po  -  P  l ) / ( po  +  pi) 

k2i  *  (p2  -  P i )/( p2  +  Pl) 

k-»o*  (p- i -  pfl)/(p.1+  Po) 


(4.21) 

(4.22) 

(4.23) 


and  doing  the  necessary  algebra,  „e  arrive  a.  ,h.  desired  solution 
for  V, .  the  potential  in  the  water  layer; 


Vf 


S' 


(FoC2+  1)  + 

'  r rr 


C  j  (  f  g  +  1 ) 


cTT 


Jo(Ar)dA 


(4.24) 


where 

-  k°*  -  Cl  -  c > ( C !  -  k0l) 

c  2  (  c  X  k  0  1  -  TJ  +  C  3  (  C  1  -  1  ) 

Ci  -  k2ie"2Ah 

C2  =  e2Ab 

C s  *  k_10  e~2X! 


(4.25) 

(4.26) 

(4.27) 

(4.28) 
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For  the  case  of  practical  Interest,  where  medium  2  is  air,  the 
source  and  measuring  elect,  ides  are  both  on  the  sea  floor, b  »  o,  l  »  0, 
and 


»  Ci  ■  e 

Fo  ■ 


-2Ah 


,  and  C2  *  1.  Then,  letting 


(k..  -  •‘2*h)(l  +  C,) 


(«'Wh  k,,  -  1)  +  C,(e'm  -  1) 


leads  to 


■£ 


Vi  »  (F# '  ’  1)  J0  (Ar)dA 


For  the  simplest  case,  where  £-►«  and  there  Is  only  a 
single  sediment  layer,  C  j-**  0  ,  and 

'k..  - 


h 


■0  "  Vo/z-o  r  +  ^ 


7™*  id .  1 


(*♦.29) 

(*♦.30) 


(Ar)dA  (4.31) 


For  the  case  where(r/l^0  (the  water  Is  very  deep,  or  the  array  very 
sma II), 


do  (°° 

*  q«Jo-k0 1  Jg(Ar)dA  , 


<»«  r  2p»  I 

r  LP«  +  pj  (4.32) 

If  (r/h)-o(the  array  is  very  large,  or  the  water  vary  shallow), 

<i«  f-rk.i  -  n 

V“*‘F"  +  qa)oir,i  +  1  J  JrCAr)dA  . 

or 

/  Zq,  Ip, 

Vr/ h—  *  —  '  (4.3?) 

which  is  the  corract  value  for  a  point  source  over  a  half  space. 


vL~, 
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A-2-4  Deep  Water,  Slngle-Leyer  Approximation 

In  many  actual  situation,.  ,h.  water  depth  ,„d  th.  thickness 
Of  the  first  layer  of  sediments  will  be  much  greater  than  the  electrode 
spacing  a.  For  this  case,  a  useful  formula  relating  the  formation 
factor,  p  sediment/  p  water,  to  the  measured  parameters  a,  AV/I,  and 
Pw  may  be  derived  from  equation  (4.32)  as  shown  below. 

From  equation  (4.32),  with  h  and  I  In  figure  4-1  much  greater  than  a 


where  pw  Is  the  resistivity  of  the  water,  and  p.  Is  the  resistivity  of 
the  sediment.  For  the  Wenner  array, 


AV  -  2(Va  -  V,a) 
or 


solving  for  Pj./pw, 


Pj/p^ 


1 


2ira 


17 

T 


This  is  the  formula  used  In  plotting  figure  2-20. 


(*♦.34) 


(*♦.35) 


a 
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A-2-5  Computer  Solution 

Equations  (4.24),  (4.30)  and  (4.31)  cannot  be  expressed  in  a 
closed  form  suitable  for  calculation  by  hand.  It  was  decided,  there¬ 
fore,  to  set  up  a  computer  program  to  solve  these  equations  by  numerical 
Integration.  Only  the  programming  for  equation  (4.31)  will  be  described 
In  detail  here.  Programs  for  the  more  general  cases  are  more  complex, 
but  similar  in  general  approach. 

The  integral  part  of  equation  (4.31)  is  In  the  form  of  a  Bessel 
function  multiplied  by  the  modulating  function  F,  where 


F 


k0l  - 

'pm 


e-2Ah 
k#  i  - 


I 


(4.36) 


V a  ■  ~  +  q0j  F(k0i.  A.  h)J,(Ar)dA  (M 

These  two  functions  are  shown  graphically  In  figure  4-2. 

J0{Ar),  approaches  zero  and  F  approaches  -k0i,  as  Ar  approaches 
Infinity. 

The  program  finds  Afl ,  the  value  of  A  at  which  F  Is  equal  to 
-k„i  wit!  i  some  small  specified  error,  it  then  evaluates  the 
Integral  portion  of  (4.31)  numerically  by  using  the  equation 


FJ0 (Ar)dA 


SA-Aa 

FJ« (Ar)dA 
A-o 


(4.38) 


The  program  used  to  calculate  the  values  of  V  for  figure  2-15  is 
included  as  section  4.2.6. 
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Ot 


F  *  -i0|  el  a  '  \a 


Not  to  scale 


Figure  4-2.  Modulated  Beeeel  function. 
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A-3  Relationships  for  Viscoelastic  Model 


A- 3“  1  Notation  Summary 


Symbol  Description 

as  Shear-wave  attenuation  coefficient 

aP  Compress  Iona  I -wave  attenuation  coefficient 

E  Young's  modulus 

Kl  Bu,k  modulus  representing  wave  energy  damping 

KR_j  8u,k  modulus  representing  elastic  response 


Shejr-wave  dissipation  function 
Compress  Iona  I -wave  dissipation  function 
Shear-wave  velocity 
Compress  Iona  I -wave  velocity 
Bulk  density 

Lame  constant  representing  wave  energy  damping 
Lamfc  constant  representing  elastic  response 
Rigidity  representing  wave  energy  damping 
Rigidity  representing  elastic  response 
Poisson's  ratio 
Angular  frequency 


A~3"2  Measured  Parameters 


4 


Compress ional -wave 
diss ipation 
function 


Shear-wave 

velocity 


i 


A  -  3  -  3  Sediment  Properties 


Property 


shear-wave 

dissipation 

function 


Symbol  Method  of  Calculation 


-1  2was 


compress  Iona  1- 
wave 

dissipation 

function 


-i 


_  i  2wa 


!/W  • 


*P‘ 


rigidity 

representing 

elastic 

response 


yR  -  0Vs2 


*  <t, 


-2 


+  2 


2(1  +  ds"2) 


rigidity 
representing 
wave  energy 
damping 


bulk  modulus 
representing 
elastic 
response 


‘W1 


- 


V* 


+  v2  +  1 

2(1  +  T£-“ST 


3  PR 


bulk  modulus 
representing 
wave  energy 
damping 

Lame  constant 
representing 
elastic  response 


■  V,(kk  ♦W  - 


i  ui 


kr  -  J  »R  ! 
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A  -  3  ~  3  Sediment  Properties  (cont'd) 


Property 


Symbol 


Lame  constant 
representing  1 
wave  energy 
damping 


Poisson's  Ratio  a 


Young's  modulus  E 


Method  of  Calculation 


A-4  Derivations  of  Rayleigh  Reflection  Coefficient 
Equations  and  Bottom  Loss 

The  intensity  of  an  echo  which  is  specularly  reflected  in  the  simple 
model  of  the  sea  floor  described  previously  is  given  by: 

>R  "  »S  '  *  *  -{J 


ILL,  -a2D 

tm2  e 


(Eq.  I) 


whe  re : 


K 

i 


*  echo  intensi ty 

=  source  intensity 

=  fractional  loss  of  intensity  at  the  sea  floor 

*  transmission  loss  due  to  spherical  spreading 


^a2D  transmission  loss  due  to  dissipative  attenuation 


furthermore; 

where: 

PR 

PS 

Pi 

ci 


of  sound  In  seawater 


Pr2 

—  (Eq.  2,  and,  I$ 


!L 

PlCl 


(Eq.  3) 


=  pressure  of  the  echo  pulse 

*  pressure  of  the  source  pulse 

*  density  of  seawater 
=  velocity  of  seawater 


Eq .  I  becomes : 


PR  P  2 

— —  -  -5 —  .  K  .  1  -  ot  2  D 

Picj  p,Cl  *  TJd) 7  e 


(Eq.  4) 


assuming  P g  is  pressure  estimated  at  D  =  I. 

Taking  the  integral  with  respect  to  time  of  both  sides  of  Eq.  4  gives: 

TTTrS  P«dt  "  pTcT  K  TOT2'  *  a2D£  ps2dt  to.  5) 


where: 


T  =  echo  arrival  time 
t  =  time  duration  of  transmitted  pulse 
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sine  wave, 


Since  the  source  generates  a 


P$2dt  -  P$2  (RMS)  •  t 


(Eq.  6) 


where  the  root  mean  square  (RMS)  of  the  pressure  is  equal  to  the  peak 


pressure  divided  by  the  square  root  of  two,  Eq.  5  becomes: 

-T+t 

p4;  \  Vdt  '  iT^T  K  Tm-  •  Ps!0>«).t  (Eq.  7) 

It  is  seen  that  K  may  be  defined  on  a  pert  pressure  basis  through  Eq.  4 


thusly; 


K  -  p‘“T  2  -a2D 

PS  • 


(Eq.  8) 


,f  PINC  represents  pressure  of  the  wave  front  incident  on  the  seafloor  and 
PREF^ep resents  the  pressure  of  the  wave  front  reflected  from  the  seafloor, 
then,  K,  the  fractional  loss  of  intensity  at  the  seafloor  is  given  by; 


.  piREF /pimc 

P i c i  /  Pici  \P INC  / 


(Eq.  9) 


For  the  model  and  incident  plane  waves,  the  ratio  of  the  pressure  of  the 
wave  front  incident  on  the  seafloor  to  the  pressure  of  the  wave  front 
reflected  from  the  seafloor  is  given  by  the  Rayleigh  Reflection  Coefficient 

(17)  (18)  thusly;  _ — 


_  VcVc,2  -  Sln*6_ 

Pi  yf\  -  S  I  n29 _ _ 

IT  v/r  e  2  /  C  2  2  -  Sln'9 


where : 


INC  +  Vci2/c2 

Pi  \}\  -  S  l  nz0 


.  10) 


d,  =  the  angle  of  incidence  (for  normal  incidence  *  -  0; 


which  was  the  case  for  this  investigation) 

The  Rayleigh  Equation  (Eq.10)  is  simplified  for  the  case  in  hand  thusly; 


P2C2  "  P  1  C  1 
P2C2  +  P  1  C  1 


(Eq.  11) 
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where  the  expression  pc 


is  called  the  acoustic  impedance  of  the 


medium  and  is  designated  as  Z . 
Z2  -  Zi 
R  “  Z2  +  Zj 


Eq.  11  becomes  by  substitution: 

(Eq.  12) 


where 


z2  =  acoustic  impedance  of  seawater 

Z1  =  acoustic  impedance  of  the  sediment 

The  relationship  between  the  fractional  loss  of  intensity  at  the  seafloor, 
K,  and  the  Rayleigh  Reflection  Coefficient,  R,  is  shown  by  Eq.  9  to  be: 

K  77  r2  (Eq.  13) 

Since  it  is  convenient  to  express  hydroacoustic  measurements  in  decibel 
form,  taking  ten  times  the  logarithm  of  both  sides  of  Eq.  8  gives; 

10  log  K  =-20  log  Pg  ♦  20  log  PR  ♦  20  log  2D  ♦  10  a  2D  log  e 


which  may  be  expressed  as 


(Eq.  14) 


BL  =  SPL  "  EPL  '  TLS  “  TLA 

where : 


(Eq.  15) 


BL  =  Bottom  loss  (in  dB)  and  is  equal  to  -  10  log  K 

SPL  =  Pressure  level  of  the  source  pulse  (in  dB//l  dyne/cm2) 
and  is  equal  to  20  log  Pg 

EPL  *  Pressure  level  of  the  echo  pulse  (in  dB//l  dyne/cm2)  and  is 
equal  to  20  log  PR 

TLg  =  spreading  loss  (in  dB)  and  is  equal  to  20  log  2D 
TLa  *  attenuation  loss  (in  dB)  and  is  equal  to  10ct2D  log  e 
Finally  taking  Eq.  13  and  stating  it  in  dB  gives: 

10  log  K  =  20  log  R  (Eq.  16) 
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which 


where; 


may  be  expressed  as: 

BL  =  20  log  R 

(Eq.  17) 


BL  =  Bottom  loss  in  d B 
R  =  Raylei«h  Reflection  Coefficient 
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A- 5  Core  Site  Locations 


266.98 


B  1 

37 

1724 

1  202.13  261.29 

2  — 

c  1 

32 

1724 

193.05 

254.67 

2 

32 

1748 

192.62 

254.96 

D  1 

30 

1830 

186.02 

251.12 

2 

31 

1900 

185.81 

251.01 

E  1 

26 

1945 

177.83 

245.56 

2 

30 

2124 

177.79 

244.55 

F  1 

27 

1706 

176.46 

259.04 

2 

27 

1720 

176.40 

259.04 

G  1 

28 

1811 

183.21 

263.04 

2 

28 

1821 

183.24 

263.04 

H  1 

32 

1759 

191.39 

268.02 

2 

32 

1816 

191.47 

268.02 

1  1 

31 

1856 

198.90 

273.08 

2 

30 

1918 

198.97 

273.07 

J  1 

36 

1802 

205.70 

276.42 

2 

37 

1818 

205.97 

276.53 

K  1 

42 

1900 

214.36 

282.31 

1 439279. 424  518634.575 
1439274.959  518635.635 
1440507.802  518442.381 


14/ 


1441868.575 

1441927.754 
1442911.054 
1442949.558 
1444141.465 
1444159.179 
1444304.657 
1444309.120 
1443305.041 
1443300.518 
1442097.923 
1442087.389 
1440994.416 
1440982.243 
1439982.869 
1439942.586 
1438707.668 


Bucket  Sample 


518126.137 
518222.873 
518143.581 
518145.261 
517888.568 
517710.842 
520229.775 
520230.355 
520380.028 
520377.636 
520561.401 
520557.192 
520792.279 
520785.110 
520761.861 
520756.106 
520986.803 
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Raydist 

Core  Sample  dw  Time  Coordinates _  Lambert  Coordinates 

Slte  No-  (ft)  (GMT)  Red  Green  x 


L 

1 

39 

1657 

215.06 

299.24 

1438892.253 

524055.829 

2 

39 

1711 

215.02 

299.19 

1438900.723 

524052.237 

M 

1 

34 

1751 

206.07 

293-92 

1440217.857 

523912.057 

2 

34 

1806 

206.01 

294.46 

1440241.186 

524020.239 

N 

1 

34 

1736 

199.09 

289.66 

1441240. 192 

523752.573 

2 

34 

1750 

199.04 

289.65 

1441251.291 

523756.840 

0 

1 

27 

1833 

191.11 

284.07 

1442385.654 

5234)5.116 

2 

27 

1846 

191.08 

284.04 

1442389.989 

523414.000 

P 

1 

26 

2029 

182.87 

280.05 

1443612.313 

523356.853 

2 

28 

2047 

182.70 

280.03 

1443639.124 

523366.384 
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Vol .  Wet.  Ga  Ga  Water.  Void  'Poros.  Gravel [Coarse' Medium  Fine  Silt  TTay  She 1 1 

Sample  Smpl .  wet  dry  Cont.  Ratio  Sand  Sand  (Silt  &  Clay 


